
by 

S. Neeraj 



DEPARTMENT OF MECHANICAL ENGINEERING 

ll\U3IAN INSTITUTE OF TECHNOLOGY KANPUR 

April, 1999 


Computer Aided Synthesis, Analysis and 
Animation of Planar Six-Bar Linkages 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 

Master of Technology 


by 


S. Neeraj 



DEPARTMENT OF MECHANICAL ENGINEERING 
INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

April, 1999 





] MAY iqqq;- 

i.iBH 
niHiS 




ik 1Z8072 


-l'V\ 

N f j:..- - 



To the memory of my beloved Nanaji 

Late. Shri V. Rajagopala Iyer 

who was known as Doctor of Machines during his Textile Mill days 


To the process of 

Creativity 

which brings something into new being 



CERTIFICATE 





It is certified that the work contained in the thesis entitled Computer Aided Syn- 
thesis, Analysis and Animation of Planar Six-Bar Linkages by S. Neeraj, has 
been carried out under my supervision and this work has not been submitted elsewhere for 
a degree. 



Dr. A. K. Mallik 
Thesis Supervisor 

Department of Mechanical Engineering 
I. I. T. Kanpur 


Acknowledgements 

I owe a deep debt of gratitude, respect and love to Mallik Sir, whose guidance and encour- 
agement made me accomplish this work successfully. He is more than a Professor to me 
- He is a Teacher who taught me to understand knowledge in its correct perspective. I 
enjoyed the care and affection he showered on me. During the entire work he allowed me to 
work at my leisure. He never set any deadlines for me and I never failed to use this undue 
privilege. I learnt many things from him - Dynamics, Mechanisms and a positive attitude 
towards life. The only thing he never taught me was ’Computers’. I pray to God for his 
health and hope for a continuing relationship with him. Working with him is a memorable 
experience for me. 

I am grateful to Dr. Sudipto Mukherjee for his guidance in the programming part of the 
thesis. I take this opportunity to appreciate Mr. Niteesh, Gaurav, Ishan and Vishesh for 
the wonderful work they have done in the project. They always helped me whenever I got 
stuck in the programming. 

I would like to acknowledge Anjanda for his help and suggestions in ’Optimizing Tech- 
niques’ in the thesis. Deep thanks goes to Santosh Kulloo without whose great help I 
would never have been able to write my report. 

Warm appreciation and regards are for Goutamda, who always bore my numerous queries 
and who was most of the time the first person to whom I would show my results. In 
addition I would like to thank all the members of the Vibrations and Dynamics Lab. I vdll 
always remember the rides on the cycle with Yogesh during which we discussed whatnot. 

I acknowledge my room partner Amreesh who was always there to encourage me. Appre- 
ciation goes to Asit and Sandeep with whom I spent a major part of my time. Thanks 
to Amit and Dhish for their ever supportive attitude. I thank Anya, Ganesh, Santosh, 
Subodh, Suma, Yogya and the whole Ghati Mandal for their affection towards me. I will 
always cherish the lighter moments of life I shared with Anjanda. Friendship with Arvind, 
Chandu, Jeyakumar, L.Shree, Santosh Kalva, Lt. Sanjay, Venkat, Shastriji is just a begin- 



ning. 

I will always remember my cute Samtron terminal. 

I wish to express my unbounded love for my Naniji, Ma, Daddy and Bhai. 
Stay at I.I.T. was a very beautiful and important phase of my life. 

Neeraj 


iii 



Abstract 


Many software packages are available for design of mechanisms. Using computers, a de- 
signer can solve complex and practical problems rapidly and at the same time can also test 
a vast spectrum of alternatives. 

The main objective of the present thesis, as a part of DST sponsored project, is to 
develop a methodology for computer aided design of Planar Six-Bar Linkages. The software 
package developed provides all the facilities of analysis and animation as well. 

The use of various optimization techniques in synthesis problems results in satisfactory 
solutions. The Least Square Error Minimizing Technique has been used which does not 
involve the time consuming iterations. Non-linearity in the design equations is resolved 
using Linear Superposition Technique. Examples are included to show the practical ap- 
plication of the synthesis methods. Computer programs have been developed using C-l— r 
(Borland 4.5 with graphics facilities). 
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Chapter 1 


Introduction 

1.1 Introduction 

A mechanism is a device that serves the purpose of transferring motion and/or force from 
a source to an output member. Mechanisms have been devised and used since the dawn 
of human civilization. Prom the modest beginning with wheel and pulley (on an axle), the 
lever and the gear, through the sewing machine and the typewriter, up to the most modern 
aircrafts and ships, the development has been tremendous. 

The early designs were totally based on imagination and intuition. The Industrial 
Revolution provided great impetus to mechanism design. Steam engines and weaving 
looms came into existence. However, it was not until the end of 19th century that a more 
systematic study of mechanisms started with ‘Theorestische Kinematic” by Reauleaux in 
1875. Since then the mechanism design has taken the form of a well established subject. 
Later on, various machine tools, textile and other heavy machinery such as bulldozers, 
lift-jacks and construction machines were also developed with the aid of ever increasing 
research in theory and operations of mechanisms. The mechanism design did not remain 
limited to the heavy machine industries only. The fields of photographic instruments, 
clocks, artificial limbs and mechanical toys also saw newer and improved designs at the 
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same time. 


The 1950s saw the introduction of digital computers. The advent of computers in the 
field of engineering, revolutionized the mechanism design. The speed and computational 
capabilities of computers rendered complex design calculations feasible and economical, 
thus giving birth to Computer Aided Mechanism Design (CAMD). 

The graphical techniques of mechanism analysis and synthesis as well as animation 
were brought under the ambit of CAMD as soon as graphics packages became available. 
The arduous task of repetitive designs to choose the better one has been rendered simple 
with CAMD. Also newer synthesis techniques, using various optimization methods, were 
developed to solve complicated problems. 

During the last 30 years, many CAMD packages have been developed. The earliest pack- 
ages were KINSYN and LINCAGES. After these came IMP, DRAM, ADAMS, MCADA, 
DADS and many others. All these packages have been quite popular in academic institu- 
tions as well as in the industry. The reason for their popularity can be attributed to the 
fact that CAMD enables the designer to use his creativity to full extent, which lifts the 
limitation on imagination leading to better and useful designs. Writing about KINSYN, 
Kaufman commented 

"In order to innovate freely in the field of kinematics, one should be unen- 
cumbered by time - consuming and distracting analytical procedures. Ideally 
a designer should be able to envision a new device and "presto"! Instantly a 
machine shop model materializes in his hands, one that he can handle, massage, 
rock back and forth, wiggle and contemplate; a model whose dimensions are 
plastic, so that he can push and pull it into new shapes, modeling it until its 
behavior is just right; a model whose topology is flexible, so that he can study 
the effects of changing a slider to a hinge, or a cam to a linkage." (Kaufman, 

R. E., Mechanism and Machine Theory, 1973, Vol. 8, p. 470) 

The objective of the present thesis is to prepare a software package, as a part of DST 
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sponsored project, for dimensional synthesis, analysis and animation of Planar Six-Bar 
Linkages. 

A mechanism which has got all the joints as lower pairs is commonly called a linkage. 
In a planar linkage, all the links move in planes that remain parallel to each other, and the 
axes of all revolute pairs also remain parallel to each other. 

Dubey [1] has prepared a software package under the same project which focuses on 
planar four-bar linkages. That package provides the tools to carry out the dimensional 
synthesis, analysis and animation of planar four-bar linkages. The synthesis is based on 
two approaches, viz., the precision-point technique and the optimization technique. The 
user has got the choice of using any of the two approaches. 

1.2 Objectives and Scope of the Present Work 

The planar six-bar linkage is a single degree of freedom linkage. It is used in situations 
where a planar four-bar linkage fails to satisfy the kinematic requirements. 

The primary objective of the present thesis is to develop a computer aided design 
methodology for planar six-bar linkages for the purpose of function generation, path gen- 
eration and motion generation. The analysis and the animation of the linkages are also 
made available to check the functioning of the linkages so designed. 

In Chapter 2, different types of planar six-bar linkages are considered. The different 
linkages that can be obtained from planar six-bar chain by fixing one link at a time, are 
discussed. Some examples of typical applications are mentioned. 

A number of optimization methods are available for linkage synthesis. The present work 
uses the Least Square Error Minimizing Optimization Technique, more popularly 
known as Least Square Method (LSM). Chapter 3, starts with a brief discussion on the 
available methods of linkage synthesis. It explains the LSM in detail. The planar six-bar 
linkage is a multi-loop linkage. The dimensions of the optimum linkage are determined 
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by minimizing the error in the loop-closure equation of each loop of the linkage. The 
multi-loop linkage gives rise to nonlinear equations which are formed as linear equations 
using the Linear Superposition Technique. The linear superposition technique is also 
explained in Chapter 3. 

The mathematical background of the synthesis procedure is outlined in Chapter 4. 
The optimum synthesis of Stephenson-I type linkage for Function Generation and Path 
Generation and optimum synthesis of Watt-II type linkage for Motion Generation are 
detailed. 

In Chapter 5, the optimization of transmission angle of planar four-bar linkage for 
Function Generation and Path Generation is given. This is an extension of the work done 
by Dubey [1]. 

The main conclusions of the present work are summarized in Chapter 6. Some directions 
for future work are also mentioned therein. 

Appendix A, gives the description of the tools of analysis and synthesis that the present 
software package provides. 

To meet the computational and animation requirements, C-l-f (Borland 4.5 with graph- 
ics facility) on DOS platform was selected as the programming language. 

The development of the package and hence the project has been done with the aim 
that the computer should accept performance specifications of the linkages in the user’s 
terms. The user should be able to test alternatives and carry out the detailed analysis 
processes without worrying about the calculations involved. At the same time, it should 
not be expected that for any given data, not in compliance with equations involved, the 
package will produce a satisfactory result. The package should be useful in both academic 
institutions and industry. 
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Chapter 2 


Planar Six-Bar Linkages 

2.1 Introduction 

The four-bar linkage has the least possible number of links of all the planar linkages. Over 
the last 100 years, various graphical and analytical techniques have been developed for its 
analysis and dimensional synthesis. The literature is very rich with numerous uses and 
examples of four-bar linkage. Still, sometimes it happens that a four-bar linkage does 
not provide the required level of performance needed for a particular application. In some 
applications a four-bar linkage fails completely. This is where the designer resorts to planar 
six-bar linkages. 

A six-bar linkage has more complex motion relationships among its links as compared 
to a four-bar linkage. Such complex relationships, either can not be fulfilled, or can only 
be partially fulfilled, by a four-bar linkage. Sometimes the four-bar linkage may satisfy the 
required kinematic objectives, but then its transmission quality may turn out to be too 
poor. 
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2.2 Classification of Planar Six-Bar Linkages 

From Griibler’s criterion, it can be easily shown that a planar six-bar linkage must have 
four binary and two ternary links. Two different arrangements are possible. In one arrange- 
ment (Figure 2.1), two ternary links are connected by a binary link, which is known as 
Stephenson’s chain. In the other arrangement (Figure 2.2), two ternary links are connected 
to each other, which is known as Watt’s chain. 


E A 



C 

Figure 2.1: Stephenson’s chain 
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By kinematic inversion, three different types of linkages are obtained from a Stephen- 
son’s chain. Stephenson-I is obtained by fixing one of the ternary links (5 or 6 in Figure 
2.1). Stephenson- II is obtained by fixing one of the binary links that are connected to each 
other (3 or 4 in Figure 2.1). Stephenson-III is obtained by fixing one of the binary links 
that are connected to the two ternary links at both ends ( 1 or 2 in Figure 2.1). 

Two different types of linkages are obtained by kinematic inversions of Watt’s chain. 
Watt-I is obtained if one of the ternary links is fixed (5 or 6 in Figure 2.2). If any of the 
four binary links is fixed (1, 2, 3 or 4 in Figure 2.2) Watt-II is obtained. 

2.3 Applications of Planar Six-Bar Linkages 

Some typical examples of planar six-bar linkages are covered in this section, which bring 
out their distinct applications. 

2.3.1 Stephenson’s Linkage 

This linkage is generally used for the following purposes; 

1. For output link oscillation with double or triple advance and return during one rev- 
olution of the input link. Figure 2.3 shows that output link ED oscillates twice 
during one revolution of the drive crank AB (Hain [2], pp. 383-5). The coupler curve 
generated by the coupler point C is indicated by a in the figure. 

2. For driving the two rocker links of a four-bar linkage, with a continuous crank rota- 
tion. Figure 2.4 shows that the revolving crank ED, via link DC, drives the double 
rocker linkage ABGF (Hain [2], pp. 386-7). 

3. For large angle of output oscillation. In Figure 2.5, crank ED makes output link AB 
to oscillate through 180®(Hain [2], pp. 390-1). 
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4. For producing dwell, where the output link remains stationary (at least approxi- 
mately) during continuous motion of the input link. Figure 2.6 shows that link ED 
has a continuous rotary motion with three dwells per revolution, with link AB as the 
driving crank (Hain [2], pp. 391-2). Figure 2.7 shows a linkage in which the output 
link ED dwells when the coupler point moves through the portion of the coupler 
curve having almost constant radius of curvature (Tao [3], p. 80). 

5. For producing a unidirectional output-link oscillation of more than 360° with a contin- 
uous rotary input (Shukla and Mallik [15]). It is seen from Figure 2.8, that clockwise 
rotation of crank BA produces oscillation of link CD of 14.5° in clockwise direction 
and 374.5° in counterclockwise direction. 

2.3.2 Watt’s Linkage 

This linkage is used for the following purposes: 

1. For obtaining a large oscillation angle of the output link from a small input link angle 
or vice versa. Figure 2.9 shows that the difference in lengths of the two arms of the 
link CGD will amplify or reduce the curvilinear displacement of the input GD of the 
second four-bar linkage GDEF thus directly varying the output angle Bqut (Tao [3], 
pp. 24-5). 

2. For obtaining a more suitable practical design. As can be seen from Figure 2.9, if 
the given hinges A and F are too far apart for a single, sturdy four-bar linkage, a 
six-bar linkage can be designed by suitably choosing the intermediate fixed pivot G 
(Tao [3], pp. 24-6). 
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2.4 Other Examples and Applications 


Many more examples and applications of planar six-bar linkage using its complicated mo- 
tion transfer characteristics are given by Main [2], Tao [3] and Hunt [4]. Erdman and 
Sandor [5] (Vol. I, pp. ,11-5) gives some typical examples of six- bar linkage as a tape 
carrier linkage in cassette player and as a prosthesis knee mechanism. Mallik et al [6] gives 
the examples of parallel-jaw pliers (pp. 364-7), fork lifter without a vertical guide (pp. 
254-8), and Hart’s straight line mechanism (pp. 372-3). These linkages are obtained from 
Watt’s chain. 

The purpose of inclusion of these examples in this chapter is to show the advantages 
of six-bar linkages over four-bar linkages. In all the cases considered, the four-bar linkage 
would have failed to perform satisfactorily. Yet, it is worth mentioning that most of the 
above examples have the origin of their design in four-bar planar linkages only. The details 
of the design and synthesis can be obtained from the references. 

The planar six-bar linkage has been a topic of research since 1950s; still a long way is 
to be covered to understand it completely. 



Figure 2.3: One revolution of drive crank AB produces two oscillations of output link ED. 
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Chapter 3 


The Least Square Error Method and 
the Linear Superposition Technique 

3.1 Methods of Linkage Synthesis 

Research in linkage synthesis is being carried out on two broad fronts - (a) Precision point 
synthesis, and (b) Application of mathematical programming techniques. 

In the precision point approach, a linkage is obtained which can generate a function, 
path or motion with exact accuracy at a limited number of discrete points. But at re- 
maining configurations there exists some difference between the desired and generated 
characteristics which is called the structural error. The precision points are selected so as 
to minimize the maximum structural error. The number of precision points is limited by 
the number of design parameters. Mallik et al [6] and Hartenberg and Denavit [7] can be 
referred to for detailed procedure of precision point synthesis. McLarnan [8] has dealt with 
precision point synthesis of planar six-bar linkages for function generation. 

The application of mathematical programming techniques to linkage synthesis provides 
a method in which the structural error, which may not be zero anywhere, is minimized 
in an overall sense in the entire range of motion of the linkage. Any number of design points^ 


^The points should not be called precision points. 



are chosen and some error quantity or function, suitably defined, using these chosen points 
is minimized. The error quantity or function is constructed using the difference between 
the actual and desired linkage characteristics. The number of design points is not limited 
by the number of available design parameters. The optimization technique to be used for 
minimizing the error can be Least Square Error method (LSM), Galerkin’s method, 
Cauchy’s steepest descent method or Random search method. 

Of all these methods the least square error method has some advantages. The problem 
formulation is very straight forward and the solution procedure is not based on time con- 
suming iterations. The nonlinear objective equations can be linearized using the Linear 
Superposition Technique. The Least Square Error Method and the Linear Superposi- 
tion Technique are explained in detail in the sections to follow. 

Fox and Gupta [9] provides a review and summary of the relationship between opti- 
mization methods (mathematical programming) and mechanism synthesis. It also gives a 
bibliography on the subject. 

3.2 Least Square Error Method (LSM) 

The method is one of the most popular methods of numerical techniques. In this method, 
the sum of squares of deviations of actual values from the desired values at some number 
of chosen points is minimized. 

The method can be clearly explained by taking an example. .A. simple linear equation 
of two variables x and y is considered. Suppose the equation is given as 

ax + by + k = 0 (3.1) 

and is to be fit through a given set of data points {xi,yi),{x 2 ,y 2 ): —,{xN^yN)- The 'un- 
knowns are a and b; k is some known constant. N can have any integral value > 2. Ideally, 
values of a and b should be such that at each given point, the left hand side of the equation 


14 



should give a value of zero, but actually at almost all points it will give a non zero value, 
say e, which is called the error at that point. For the error to be equally distributed the 
sum of errors at each design point, i.e., is minimized with respect to a 

and b. Writing in mathematical terms, 

-^{el + el + ... + e%)=0 

^ (Z + by + kf) = 0 

N N N 

+ + = ^ (3.2) 

similarly, 

— {e\ + 63 + ... + e^] = 0 

^ (E = 0 

N N N 

=^aY,xy + bY,y^ + k'Ey = ^ (3.3) 

Equations (3.2) and (3.3) form a set of simultaneous linear equations which can be 
solved by various methods like Gauss Elimination method etc., to give the values of a and 
b. 


3.3 The Linear Superposition Technique 

The equation considered in the previous section, i.e., Equation (3.1) is a linear equation. 
If the equation has a nonlinear term, then LSM can not be applied directly. For example, 
an equation with a nonlinear term is taken as 

ax-{-by + abki + fca = 0 (3.4) 





and the same N data points are to be used. ki and k 2 are some constants.The term ab 

brings nonlinearity in the equation. The linear superposition theory, as put forward by 

Bagci and Lee [10], partitions the above nonlinear equation into two linear equations. This 

is shown in the following steps. The term ab is represented by a so-called nonlinearity 

coefficient A , i.e. 

ab = X. (3.5) 

Further, it is defined that 

a = pa + Xqa 

(3.6) 

b = pb + Xqb 

(3.7) 

resulting in 

{Pa + A^a) X -t- {pb + Xqb) y -h AA;i + k 2 = 0 


or 

{PaX + Pby + k 2 )+X{qa + qby + ki) =0. 

The above equation is partitioned into two linear equations as, 

(3,8) 

{PaX -1- PbV + /C 2 ) = 0 

(3.9) 

and 

{qaX + qby + ki)=0. 

(3.10) 


Both the resulting equations are solved using the LSM, giving values of (pa, Pb) and 
{Qa, qb)- The value of A is obtained from Equations (3.5), (3.6) and (3.7) as, 

{Pa + A^a) X (pfc + Xqi) = A 
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or 

X'^Qaqb + A [pbQa + PaQb “ 1 ) + PaPb = 0 . (3.11) 

The above quadratic equation, also known as the compatibility equation, is solved to 
obtain the value(s) of A. No solution exists if the discriminant of Equation (3.11) turns out 
to be negative; a single value exists if it is equal to zero and two distinct solutions if it is 
greater than zero. Once A is known, values of the original unknowns a and b are obtained 
from Equations (3.6)and (3.7). 

3.4 Use in Linkage Synthesis 

When the LSM is applied to Linkage Synthesis, the simple equations viz., Equations (3.1) 
and (3.4) are replaced by the loop closure equations. The unknowns, instead of a and b, 
are the dimensional parameters of the linkage. The data points are replaced by functions 
of the displacements which are to be generated by the linkage at the N design points (or 
positions). The error to be minimized is the sum of squares of error in the loop closure 
equations, at each design point. 

Depending on the geometry and the number of unknown dimensions of the linkage, the 
loop closure equations may be linear or nonlinear. The nonlinearity is removed by using 
the linear superposition technique. The LSM results in simplified formulation and non 
iterative solution of the design equations. 

The synthesis technique employing linear superposition has been applied to single- 
loop linkages for function generation [10], for path generation and rigid-body guidance (or 
motion generation) [11], and to multi-loop linkages for function generation [12] and for 
path generation and rigid-body guidance [13]. The present thesis presents the optimum 
synthesis techniques as applied to planar six-bar linkages for function generation, path 
generation and motion generation. The details of these are given in Chapter 4. 


-1 m 


Chapter 4 


Synthesis of Planar Six-Bar Linkages 

4.1 Classification of Synthesis Problems 

The present work focuses on dimensional synthesis of six-bar linkages. Dimensional syn- 
thesis involves the determination of link-lengths of a linkage so as to fulfill the prescribed 
motion characteristics. Depending on the motion characteristics to be satisfied, the syn- 
thesis problems can be classified as listed below: 

1. Function Generation: In this type of problems, the motion parameters (displacement, 
velocity, acceleration etc.) of the output and input links are to be coordinated so as 
to satisfy a prescribed functional relationship. 

2. Path Generation: Here a point on the floating link (i.e., a link not connected to the 
frame) of the linkage has to be guided along a prescribed path. The generation of 
the path may or may not be coordinated with the input motion. 

3. Motion Generation: A rigid body (i.e., one link of the linkage) is guided in a pre- 
scribed manner. Here also the guidance may or may not be coordinated with the 
input motion. 


18 


The subsequent sections explain function generation and path generation for Stephenson-I 
and motion generation for Watt-II type linkages. 

4.2 Function Generation using Stephenson-I Linkage 

Figure 4.1 shows a Stephenson-I linkage. The input angle 9 tn represents the independent 
variable x, while the output angle dour represents the desired function y = f{x). The link- 
vector diagram corresponding to Figure 4.1 is shown in Figure 4.2. The design variables are 
Ti2, L22, L31, L41, Tso, Lsb, Lea, Let, Tcand rjc^. The values of 0 /jv and 9out are represented, 
respectively, by 9i2 and ^41. The values of these angles at the N design points are evaluated 
using the prescribed function relationship, chosen scale factors (between x, y and 9i2 , 9 ^i ) 
and the respective initial values (of x, y and 9x2 ; ^4i)- Out of the ten link design parameters, 
one link-length is arbitrarily chosen to fix the overall size of the linkage. 



For the loop ABODE A, the loop closure equation is 


~^12 + ^56 + ^31 — ^66 ~ ^41 — 0 


( 4 . 1 ) 


^In all the sections of function generation, using Stephenson-I linkage, Lu should be taken as £i, L22 
as L2, L21 as £3 and L^iss £4. 
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Figure 4.2: Link-vector diagram of Stephenson-I at the i*^ design position for function 
generation. 

and for the loop ABGFA, the loop closure equation is 


~ ^6a — ^22 — 0 


(4.2) 


Elimination of angle ^31 from Equation (4.1), gives 


■^12 "I" -^56 d- L41 + Lgj — L31 + 2.Li2-L5b cos {0i2i — ^ 5 ai ~ 7 c) 


— 2.Z/i2.I'66 cos 6i2i — 2.L12.L41 cos (^I2i — ^41i) 2.L^b-Li\ cos (05oi 7c “ ^4i) 


2.L56*-1^66 COS (^5c^ 4 * 70 ) 4“ 2.Z/41 COS ^4|j 0 


(4.3) 


at the i^^ design position and elimination of $22 form equation (4.2), gives 


^12 d" -^Sa d" Lgg — L22 2.Li2.I'5o COS (^i2t — ^5ai) “ 2.Z/i2-L6a COS (0i2i ~ Vc) 


2.L5a*-f^6a COS (^5ai ^c) H * 


(4.4) 



The values of ii 2 , L^b, jc and rjc are assumed whereas those of 6 i 2 i and are deter- 
mined from X and y for the N design points. The values of Obai are assumed arbitrarily at 
the design points. The optimum dimensions of the linkage are those which satisfy Equations 
(4.3) and (4.4) at N design points in the least square sense. 

Dividing Equation (4.3) by 2 .L 12 .L 66 and rearranging, the input-output relationship at 
the i*^ design position is obtained as 


7 ^ (cos (041i) - ^ cos ( 0 i 2 i - 041i)) (- cos (05ai + 7c) + ^CCOS OS {dui - 9bai 

^12 \ Leb / Ti2 \ Leb 

Ll^ + Ll + Ll+Ll-Ll 


2 .Li 2 .Leb 


— ^4ii) 

)•&(- 

-ill 

Tee T 41 

Ti 2 Tee 


7c 


cos {Osai + 7c - ^41i) - cos {Oui) = 0 , 


or 


DiiJ,u + D 2 H 2 i + — DiD 2 pi — — 0 


(4.5) 


where, 

= fc . -02 = fe. O 3 = A = = fe-fe 

Mli = cos (041i) - cos [Bui - ^41t)> 

/^ 2 i = — cos (Bsai + 7c) + xtt ~ ^5at — 7c) j /^Si = Ij 

Pi = ^ COS (&5ai + 7c — ^41i)) = COS (^I2i) ■ 

Similarly, dividing Equation (4.4) by 2 .Li2Lba and rearranging one gets 

Tea \ , Tea Tf2 + Tg^j -I- T|^ — T22^ 

COS {O^ai Vc) Y COS {Oi2i 

^12 Z.lji2>Lf^a 


- cos {Bui - Bsai) = 0 , 


or 




(4.6) 


where, 




1 Li2^ 2 ’ 3 2.Li2-L5a 



= cos {Osai - Vc), l^ 2 i = COS {dm - Vc), fi'si = -1 
v'i = - cos {6m - d^ai) ■ 

The design parameters Di,D 2 ,D 3 and D[,D' 2 ,D'^ are defined in such a way that a 
design parameter appears only in either Equation (4.5) or Equation (4.6). In other words, 
both the equations are decoupled in terms of unknown design parameters. 

Equation (4.5) is a nonlinear equation while Equation (4.6) is a linear equation. Equa- 
tion (4.5) is partitioned into two linear equations using the linear superposition technique, 
as explained in Chapter 3. 

From Equation (4.5), 

A = A^2. (4.7) 

It is defined that, 

Di = Xqi, D 2 = P2 + Xq2 
and 

D 3 = P3 + Xq 3 . (4-8) 

Equations (4.8) and (4.5) give the following two linearized equations: 

PlPli + P2P2i + PzP'Zi -Vi = 0 (4.9) 

and 

QiPii + q2P2i + qsPzi - pi = 0 . (4.10) 

Equations (4.6), (4.9) and (4.10) are used to get the required unknowns, using the LSM 
for the N design points. Equation (4.6) gives D\,D' 2 ,IJ> 3 , Equation (4.9) gives Pi-P 2 ,Pz 
and Equation (4.10) gives qi,q 2 ,qz- The compatibility Equation (4.7) is expanded to give, 

qiq 2 ^ + {Piqz + P 2 Qi ~ 1) + P 1 P 2 = 0 (4.11) 
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which is solved for obtaining the value(s) of A. Once A is known, Z?i, D 2 , D 3 are calculated 
using Equation (4.8). From Di^D^^Dz and D'lyD^jD'z all the unknown link dimensions 
can be determined. 

A six-bar linkage has four geometric inversions. Thus, with two real values of A there are 
in total, eight geometric inversions. One of these eight is the optimum solution. Further, 
the resulting solution generates values of dsai which are different from the initially assumed 
values. These generated values can be used in place of the initial values to improve the 

V 

accuracy of the generated function. Such improvement can be employed a few times. 
Different values of the other assumed dimensions may also be tried to increase the number 
of acceptable optimum linkages. 

The error in the generated function at the i^^ design point (position) is 

e* = (04u)p - (^4u), (4.12) 

where ( 04 ii)^ is the generated value and (^ 4 ii)^ is the desired value. The efficiency of the 
generated approximation is tested by the average error magnitude 



and/or by the root-mean-square error 

= (4-14) 

4.3 Path Generation using Stephenson-I Linkage^ 

The Stephenson-I linkage at the i‘^ design position is shown in Figure 4.3. Coordinates of 
the design points along the path are defined by = Pilfi in the OXY system, and 'PiS 

^In all the sections of path generation, using Stephenson-I linkage, L12 should be taken as Li, L22 ^ 
ij2) -^3 and i^ 4 * 
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are specified in a design problem, cfj = cos 6ii + sin 5ij is the unit vector positioned by 
the angle Si, where i and j are unit vectors along X and Y respectively. The fixed link 
is positioned by = Rs^aB = Rp^ap- Dimensions of the 

linkage are defined in Figure 4.3. and (where /? = ^31 4- define the position of the 
coupler point P with respect to the coupler link. The loop-closure equation for the coupler 



Figure 4.3: Stephenson-I, as Path Generator, at the i*^ design position when generating a 
path. 

point, using the loop OABCPO is. 


+ ^12 + ~^5b 4* ^ — 0 ■ 


( 4 . 15 ) 
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The other two loop-closure equations, using ABGFA and ABC DBA, are 

L 12 + L 5a ~ B 22 — ^ (4-16) 

and 

^12 + -^'56 + ^31 ~ ~ ^41 = 0 ■ (4-17) 

Elimination of angle 13 from Equation (4.15) at the i‘^ design position gives 

^A + ^12 d" -^56 3- Pi — + 2 .R^.Li 2 cos(q;^ — 6i2i) + 2 .i 2 ^.X' 5 j cos(a;>i — dsai — 7 c) 

-2.RA.PiCOs{aA - 4) + 2.Ti2.L56COS(0i2i - Ooai - 7c) - 2.Li2.Pi cos(^i 2 i - 6 i) 

- 2 .L 55 .Pi cos(^5a,- -f 7c - <5i) = 0 . (4.18) 

Elimination of 022froni Equation (4.16) at the i*'^ design position gives 

■^12 + -^L + J?!' + — T22 + 2 .Li 2 .L 5 a cos(0i2i - ^ 5 ai) ” 2.Li2-Rf COS(0i2t — OLp) 

-|-2.Li2.i?A cos(012t — a^) — 2.L5a.i?ir COs(^5aj — O-p) + 2.L5(i.i?^ C0s(^5at — CK^l) 

—2.Rf-Ra cos(q:f — a a) = 0 . (4-19) 

Similarly, elimination of On fro Equation (4.17) gives 

•^12 d" -^56 d" Ln + R% d" P^ ~ -^41 d- 2. £12.1/56 COs( 012 t — ^Soi ~ 7c) 

-t-2.Li2-£'3i COS(012i “ ^31i) ~ 2.Li2.J?E COs(012i “ O^e) A 2.Li2.i?^ C0S(^12i — OCa) 

+2. £56.1/31 cos( 05 ai d- 7c “ ^Sli) - 2 .L 55 .RE COs( 05 ai + 7 “ <^e) + 2 .L 55 .RA COs( 05 ai H" 7 ” ^a) 

- 2 .Lzi.Re cos(^3ii - as) -f 2 .£ 3 i.jR^cos( 03 ii - o;^) - 2.Re-Ra cos(o;£; - a^i) = 0 . (4.20) 
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Equations (4.18), (4.19), and (4.20) are the design equations that must be satisfied at 
the N design points. At each design point values of 621 % specified. The 

values of ^311 to be used in Equation (4.20) are obtained after solving Equation (4.18) and 
using the separated x and y components of the loop-closure Equation (4.15). Thus, 


tan (03it + 1 \)) 


Pj sin(^i) - Ra sin (g^) - Ln sin (^121) - Ub sin {Bsqj -t- 7 c) 
Pi cos ( 5 i) - Ra cos (cxa) - Lu cos (Om) - Lsb cos (O^ai + 7c) 


(4.21) 


The values of L12 , q;a, a/?, 1/) and 7c are initially assumed. The remaining dimensions 
L 22 , Lzi, T41, Lba, T56, Lba, Leb, Ra, Re and Rp are obtained by applying the LSM. 
Equation (4.18) is rearranged as 


Difiii -i- D2fi2i + DzfJ'Zi — DlDlPi — I'i — 0 


(4.22) 


where, 


= Ra, D 2 — Lbb, D3 = R\ + Li 2 + L^b — a^, A = D 1 D 2 — RA-L^b, 
pu = 2.L12 cos {dui - a a) - 2 .Pi cos (a^ - <Ji) , 

P 2 i ~ 2.i!(i2 cos {di 2 i ^5ai 7c) ^-Pi COs(05ai -f" 7c ^i), PZi ~ 1) 

Pi = -2 COS (q:a - dzai - Ic), = ‘^■Lu.Pi cos (^I 2 i - 5i) - Pf . 
Equation (4.19) is rearranged as 


^iPli "b ^2p2i "b ^zP'Zi ^1^2Pi ® 


where, 


D[ = Ua, P>2 = '^E, D', = L\2 + LI + Rl + R\- LI 2 , A' = D;D; = Lza-RA, 

Pn — 2.Z/12 cos {0i2i ^oci) “b 2. COS (05ai CH/i), 

P 2 i = 2.L12 cos { 0 i 2 i - otp) - 2-i?yi cos(aF - a^i), Mat = 

Pi = 2cos(05ai - a^), v'i = -2.Li2.i?A COS (012i - OLj^ ■ 


(4.23) 


26 



Similarly, Equation (4.20) gives 


_ // n // tt ft ri ii tf 

DitJ-ii + i52^2i + ~ DiD2Pi — — 0 


(4.24) 


where, 

Di = 1 ^ 31 , £>2 = -^3 — -^12 + -^56 + -^31 + ~ -^ 41 ’ ^ = D-^D^ = L-^-Re, 

Pli = 2.Z/12 cos (0i2i ~ ^31i) + 2.L56 COS (^Soi + 7c “ ^31i) + 2.i?A COS (dzii ~~ 0 !a)) 

/i2i = -2.L12COS (0i2i - as) - 2.L5JCOS (^5ai + 7c - CTe) “ 2. i?A COS (()!£; - Q^), Pzi = 1> 
p" = 2 cos( 03 ii-Q:£), 


v'l = -2.L12.L56 cos (6>i2i - ^5at - 7)-2-Li2.i2.1 COS (0i2i - aA)-2.Lzh.RA cos ((95ai + 7c - a^) . 


Equations (4.22), (4.23) and (4.24) have nonlinearity coefficients as A = D 1 D 2 , X' 
= D'iD '2 and a" = I?! £> 2, respectively. These are linearized by the linear superposition 
technique as had been explained in the previous section for the function generation problem. 

The maximum number of solutions possible in the path generation problem is2x2x2x4 
= 32. Out of these one or more solutions will satisfy the problem up to the desired accuracy. 
The errors in the approximation are defined as follows to check the efficiency of the method. 
The average error is given by 


I fN ^ 

Bavg = ^ + {Vg - yd\)i^ 


and the root-mean-square error 




(4.25) 


(4.26) 
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4.4 Motion Generation using Watt-II Linkage^ 

Figure 4.4 shows a Watt-II linkage with link DE as the carrier link, coupler point P as the 
generator of the path. The rigid body, shown in the Figure 4.4, is positioned by positioning 
the body-line PQ. When guiding a rigid body, the angular dimension e controls the position 
of the rigid body-line with respect to the link-line of the generator link. If e is taken as 
an unknown dimension, it increases the number of compatibility equations to be satisfied, 
hence it is preferable to specify e. In the following developments, the value of e is taken to 
be known and is taken to be positive^ . At N design positions the value of (where 
= Pitf i in the OXY system, as defined in the previous section, if, = cos Sii + sin Si j is 
the unit vector positioned by the angle 5i.) and (anexit 

gle made by P^ at the design position) are specified as required by the design. The 
three loop-closure equations for the loops OAFPO, OBCGAO and FGDEF, respectively, 
are 

+ ^66 + - 'Pi = 0 (4.27) 

Pb + P 2 I + ^5a ~ ^6a ~ Pa — 0 ’ (4.28) 

and 

Psb ~ Psa + Psb ~ Pba + ^31 + ^41 = 0 • (4.29) 

It is seen from the Figure 4.4, that 

9^u = ^i-e (4.30) 


and 

Pi = ^ 41 * + (ISO** — l{}) 

®In all the sections of motion generation, using Watt-II linkage, Ln should be taken as Li, L 21 as 
L 31 as L 3 and Luas L 4 . 

^i.e., the value of c is specified between 0° and 360®. 
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or 


ft = ?(-£+ (180” -</.) . 


(4.31) 


The values of Lea, ci.-i, OiB-, 7 c, Vc and e are specified. The values of Osai are assumed for 
each design position in accordance with Pi and The remaining unknown dimensions are 
found from the synthesis. 

Elimination of angle O^b from Equation (4.27) at the design position gives 

R\ + + Pi — Lib + ^-Ra-o- cos(q;^ - Pi) — 2.RA.Pi cos(o;a — <^i) — 2.a.Pi cos (Pi - Si) = 0 . 

(4.32) 

Elimination of 02ifrom Equation (4.28) at the design position gives.R_ 


R% + Rta + + R\ ~ -^21 + S.Rs.Lsa COs(q;b ~ ^5ai) ~ ‘^-Rs-Lsa COs(arB — ^gai) 


-2.Rb.Ra cos (qb - oca) - 2.Lba.L^a cos ( 05 oi - 06oi) “ 2.Lba..RA COS {dbai “ 0!a) 

+2.Lea.RA cos{9eai - cxa) = 0 . (4.33) 

Similarly, elimination of ^31 from Equation (4.29) gives 

Reb + Rsa "k Rbb + Rla + -^41 “ -^31 ~ 2 . Left. Lea COs(r 7 c) + 2 .Leb.Lbb cos (Sgat — Vc — 9 bai — 7 c) 

2 .L(,b.Lba COS(^6oi V ^ 5 ai) T ‘^•L^b-Ln COS(06ai Vc ^ 41 t) 2.L6a-L5j COs(^6at ^oat 7cj 

T2.Lea.L5a COs(^eai ^5az) 2.Lea*L4i COs(^eai ^41i) 2.L55.L5a COs(7c) 

T2.L55.L41 COs(05at T 7c — ^4l) ~ 2.L5a.L4i COS(05ai — 041i) = 0 . (4.34) 

The values of 6 eai to be used in Equations (4.33) and (4.34) are obtained after solving 
Equation (4.32) and then using the separated x and y components of Equation (4.27), 


.^n 



Thus, 


tan {deai - Vc) 
Equation ( 4 . 32 ) is rearranged as 


Pi sin((5i) — Ra sin {a a) — a sin {Pi) 
Pi cos (Si) — Ra cos (a^) - a cos (Pi) ' 


( 4 . 35 ) 


+ DzfJ>zi — DiD 2 Pi — Vi = Q ( 4 . 36 ) 


where, 

Dx — Ra, D 2 = a, D 3 = Ra + A = D 1 D 2 = 

Pii = - 2 .PiCos (aA - Si), 

P 2 i = -2.Pi cos(Pi - (5i), pzi = 1, 

Pi = -2 cos (a A - Pi), Vi - -Pi . 

Equation ( 4 . 33 ) is rearranged as 


DiPii + D'2P2i + DzPsi - D-^D^Pi - = 0 


( 4 . 37 ) 


where, 

D[ = Rb, D'^ = Lza, P; = Rl + Ll + Ll + Rl- L^, X' = D[D'^ = R^.Lsa, 
p'li = -2.Ra cos (cXa - O^b) + 2.Lza cos (aA - 06ai), 

P 2 i = -2.Lea cos (Osai - Oeai) - 2.Ra cos(^5ai - a^), p'zi = 1, 

Pi = -2 cos (as - Osai): ^i = -2.L(,o.-Ra cos (Q^ai - OLa) . 

Similarly, Equation (4.34) gives 


u H // // it It ft it ll 

^2 l^ 2 i ^ ~ ^2 Pi 


Ui =0 


( 4 . 38 ) 


where, 

D'; = Ls,, d; = Uu d; = li + q, + if, + li + q, - q,, v = d';d: = 

Pli ~ 2 .Z/ 65 COS {O^ai “ ^5ai ” 7c) (06ai ^5ai 7c) ” 2.Z/5 q COS (7c)' 

P' 2 i “ 2.Z/5J^C0S {O^ai ^ 41 i) 2.L5^COS {O^ai ^ 41 i) “ 2.X5Q COs(05Qi — * ^411)7 Pzi 



p- = -2 cos {05ai + 7c - ^41i), 

U- = -2.L66.L6a COS {rjc) - 2.LQb.L<^a COS {Q&ai - Vc - dsai) + 2.L^a-L5a COS {6(,ai “ O^ai) ■ 

Equations (4.36), (4.37) and (4.38) have nonlinearity coefficients as A = D 1 D 2 , A — 
D[D '2 and A" = D'^D'^, respectively. These are linearized by the linear superposition 
technique. 

All these equations give maximum 2x2x2 = 8 sets of dimensions, therefore 4x8 = 32 
linkages are defined by geometric inversions, out of which one or more are the optimum 
solutions. The angular error in the generated motion at the i‘* design position is 


Ci — (^4ii + e)g - (6)d 


(4.39) 


where (^ 4 ii)p is the generated value. The average error magnitude is given by 



and the root-mean-square error as 



(4.40) 


(4.41) 


The errors in the generated path are defined by Equations (4.25) and (4.26). 


4.5 Numerical Examples 

Computer programs in C++ (Borland 4.5) have been developed which perform the opti- 
mum synthesis of Stephenson-I for function generation and path generation problems and 
Watt-II for motion generation problems. This section illustrates some examples giving the 
results generated by those programs. 


^19 



Example 1: Stephenson-I linkage is synthesized to generate the symmetric function 


y = x^, -1 < X < +1 . 

within the intervals® = 90° and A 04 i = 60°. Nineteen uniformly distributed design 
positions are used. Values of x and y are given in Table 4.1, where the subscript d designates 
the "desired value" while g designates the "generated value". The specified dimensions 
(Figure 4.2) are L 12 = 20 cm, Let = 40 cm, tJc = -142°, 7 c = —154°, ^12 = 

9li = 10°. After solving the design equations, two distinct values of A are obtained as 
—139.606 and —0.875. The optimum linkage® was obtained using A = —0.875. The average 
error is Cavg — 0.622° and the root-mean-square error is erms — 0.768°. Table 4.2^ gives 
the displacements of the linkage and the error in the generated function along with the 
specified values of ^i 2 t,^ 5 ai and dm. The dimensions of the linkage and the transmission 
angles pi and fj ,2 are also given. The values of all the link-lengths have been made positive 
with the corresponding changes in their respective angular displacements. 

Table 4.1 also gives the generated values of y at the design positions. Table 4.3 gives 
the values of the displacements obtained when the values of $ 5 ai generated in the first 
optimization process ((^sat)^ in Table 4.1) are used as the input for a improved solution. 
It is seen that the solution obtained does not prove to be a better solution. Hence, the 
linkage obtained in the first optimization process is taken as the optimum solution linkage. 

From Table 4.2, it is seen that the error in the generated function has a maximum value 
of —1.794°. In reference [7] (Table 10.3, p. 314), the maximum value of error generated by 
a planar four-bar linkage, for the same function and same input-output range of motions 
with five precision points, is given as 2.34°. Though the improvement obtained by the 
six-bar linkage is not much but it can be hoped that a better solution can be obtained with 

’^Relationship between x and Bn, and y and d 4 iis given at the end of Example 1. 

*The linkage obtained using A = —139.606 gave eavg = 24.870° and Crm* = 28.315°. 

^AU the angles, in the given Tables, are in degrees. All the numerical values are given up to the third 
decimal place. 




a different set of initial data. 

A linear relationship is assumed between Ax and given by 
«;2 - = ( i ‘ - 

Similarly, a linear relationship is given between Ay and A 04 i is given by 



the pages of a music book automatically (reference [5|, Problem 1.42, p. 67). The linkage 
is shown in Figure 4.6. To the upper right-hand corner of each page a magnetic strip is 
attached. The linkage is to have a crank input (possibly actuated by a foot switch). 

Eighteen design positions are used whose co-ordinates are given in Table 4.6. The 
specified dimensions® (Figure 4.3) areLia = 60 cm = 180°, = 0°,aF = -103°, Tc = 

—20° and ip — 60°. Only 16 geometric inversions out of 32 had mobility within the zone 

of motion. The linkage with the best performance is shown in Figure 4.7. The values of 

\ 

errors generated are eavg = 2.746 cm and erma — 2.382 cm. The minimum values of the 
transmission angles are = 25.011° and /zzmin = 32.436°. The errors in the generated 
coordinates are given in Table 4.6. Table 4.7 gives the dimensions and the displacements 
generated by this linkage. 



(b) 


Figure 4.6: Page turning linkage. 

Example 4- Watt-II linkage is synthesized to guide the door-board of a vehicle garage 
(reference [6], Problem 6.3-2, pp.246-7). The linkage is designed so as to have the motion 

of the door without touching the vehicle that is parked inside the garage, i 
®The dimensions of the linkage correspond to a scaled model and not to the actual linka 
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Figure 4.7: Stephenson-I linkage optimized in Example 3 to generate path for the page 
turning linkage. 

Figure 4.8. Both the desired (marked by o) and the generated curves traced by the end 
of the door-board are indicated in Figure 4.8. Sixteen design positions of the rigid-body 
(i.e., the door-board) are taken as given in Table 4.8. The specified dimensions® (Figure 
4.4) are Lsa = 15 cm, a>i = 0°, qb = 44®, 7 c = — 90°, rjc = 0®, 0 = 9® and e = 351®. 
Only four geometric inversions out of 32 had mobility within the specified zone of motion. 
The optimum solution, shown in Figure 4.9, guides the door-board with angular Cavg = 
1.713°, angular Crms = 2.110°, linear Cavg = 3.075 cm and linear e^ms = 2.995 cm. Table 
4.9 gives the values of the specified along with the generated displacements. 

It should be noted that the value of link-length is very small as compared to other 
link-lengths, suggesting that one can make the hinge C of second order where links 2, 3 
and 5 are all connected. Table 4.10 and Table 4.11 give the values of the displacements of 
all the links under such a situation. It may be mentioned that the door-board, taken as 
the coupler of a four-bar linkage could never achieve this motion. It would have required 
a much larger overhead space in order to avoid the vehicle body. 

^The dimensions of the linkage correspond to a scaled model and not to the actual linkage. 
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Figure 4.8: The available space inside the garage is determined by the envelope of the 
straight lines, representing the garage-door at various positions. 



Figure 4.9; Watt-II linkage optimized in Example 4 to guide the garage-door. 



4.6 Conclusions 


The Least Square error Method and the Linear Superposition Technique provide a powerful 
and efficient synthesis method. The number of unknown dimensions can be increased in 
each case of synthesis. This in turn increases the complexity of the equations leading to 
more than one compatibility equations. For the solution of more than one compatibility 
equations, Sylvester’s Dialytic Method of Elimination^ has to be used [11-13]. The number 
of solutions obtained also increases with the increase in the number of unknowns. Hence to 
keep the mathematics simple and the number of solutions manageable, the synthesis with 
higher number of unknowns has not been pursued in the present work. 

Different initial values of the specified variables result in different linkages, thus permit- 
ting infinite number of solutions. During the calculation, conditions may arise, such that 
any parameter gets a non-real value. In such conditions, the synthesis should be stopped at 
that place and restarted with changed initial values. However, this choice of initial values 
is very critical for obtaining a good solution. A number of trials is often necessary to arrive 
at an acceptable solution which can then be easily refined to obtain a better solution. This 
shortcoming becomes more pronounced, when a crank is needed as the input link (contin- 
uous rotation in one direction). The difficulty can be well understood by taking the case 
of function generation. Referring to Figure 4.1, suppose link 1 of the Stephenson-I is to 
be a crank for certain oscillation of the rocker link 4. The problem becomes complicated 
because the loop ABGF, is now a four-bar crank rocker, hence the rocker, i.e., link 2 can 
not cross the line of frame, i.e., the line made by AF. Moreover the loop ABGF should 
satisfy the Grashofian condition, with link 1 as the smallest link. The initial specifications 
of Tefc, L 5 („ 05 a, 7 c and r]c should be such so as to satisfy the above conditions with min- 
imum error in the generated function. The complexity further increases in the cases of 
path generation and motion generation, where there are three loop closure equations to be 
satisfied, leading many times to abandon the initial set of the specified dimensions and the 


^For details see Bocher, M., Introduction to Higher Algebra, Macmillan, 1949, pp. 195-202. 



angular values. 

Apart from the values of the error in the approximations, transmission angles as well 
as the maximum link-length ratio (i.e., compactness of the linkage) can also be taken as 
criteria for selection of the optimum 

The Least Square Error Method for loop-by-loop synthesis of planar six-bar linkages 
for function generation, path generation and motion generation, using the linear superpo- 
sition technique leads to efficient linkage designs without involving iterative solutions. The 
process can be repeated by reusing the generated value to improve the efficiency. 
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Table 4.1 Desired displacements for the function ( y = -l<x<+l) for 

which Stephenson-I linkage is optimized in Example 1. 


i 


(yih 


{yi}g 

1 

-1.000 

1.000 

1.005 

0.005 

2 

-0.889 

0.790 

0.792 

0.002 

3 

-0.778 

0.605 

0.609 

0.004 

4 

-0.667 

0.444 

0.448 

0.004 

5 

-0.556 

0.309 

0.309 

0.001 

6 

-0.444 

0.198 

0.193 

-0.005 

7 

-0.333 

0.111 

0.100 

-0.011 

8 

-0.222 

0.049 

0.035 

-0.015 

9 

-0.111 

0.012 

-0.001 

-0.013 

10 

0.000 

0.000 

-0.006 

-0.006 

11 

0.111 

0.012 

0.016 

0.004 

12 

0.222 

0.049 

0.062 

0.013 

13 

0.333 

0.111 

0.128 

0.017 

14 

0.444 

0.198 

0.212 

0.015 

15 

0.556 

0.309 

0.314 

0.006 

16 

0.667 

0.444 

0.437 

-0.008 

17 

0.778 

0.605 

0.584 

-0.021 

18 

0.889 

0.790 

0,770 

-0.020 

19 

1.000 

1.000 

1.030 

0.030 


Table 4.2 Desired and generated displacements for the Stephenson-I 
linkage in Example 1 for fxmction generation. =90", A04i =60". 


i 

^12i 

{^s<u)a 

{^ Sal)g 

{d21i)g 



(^4u)^ 


Mu 

Mu 

1 

50.000 

-170.426 - 

170.510 

55,486 

-75.940 

10.000 

9.724 

-0.276 

45.996 

85.664 

2 

55.000 

-169,012 - 

169.056 

58.998 

-73.743 

22.593 

22.451 

-0.142 

48.054 

83-806 

3 

60.000 

-168.013 - 

168.076 

62.623 

-72.146 

33.704 

33,488 

-0.216 

50.699 

74.366 

4 

65.000 

-167.452 - 

167.514 

66.338 

-70.902 

43.333 

43.110 

-0.223 

53,852 

65.988 

5 

70.000 

-167.302 - 

167.315 

70.124 

-69.821 

51.481 

51.432 

-0.049 

57.439 

58.747 

'6 

75.000 

-167.502 - 

167.434 

73.962 

-68.738 

58.148 

58.435 

0.287 

61.396 

52.827 

7 

80.000 

-167.983 - 

167.836 

77.834 

-67.501 

63.333 

63.991 

0.658 

65.670 

48.508 

8 

85.000 

-168.682 - 

168.496 

81,723 

-65.971 

SI ,031 

67.909 

0.872 

70.219 

46.120 

9 

90.000 

-169.560 - 

169.397 

85.614 

-64.047 

69.259 

70.034 

0.775 

75.011 

45.919 

10 

95.000 

-170.610 - 

170.534 

89.493 

-61.689 

70.000 

70.356 

0.356 

80.027 

47.954 

11 

100.000 

-171.857 - 

171.909 

93.344 

-58.922 

69,259 

69.027 

-0.232 

85.253 

52,051 

12 

105.000 

-173.353 - 

173.536 

97.150 

-55.808 

67.037 

66.278 

-0.759 

89.314 

57.914 

13 

110.000 

-175.173 - 

175.438 

100.894 

-52.414 

63.333 

62.318 

-1.015 

83.669 

65.268 

14 

115.000 

-177.406 - 

177.653 

104.552 

-48.795 

58.148 

57.269 

-0.879 

11 .195 

73.936 

15 

120.000 

179.863 

179.761 

108.097 

-44.991 

51.481 

51.141 

-0,340 

71.663 

83.868 

16 

125.000 

176.564 

176.716 

111,492 

-41.030 

43.333 

43.811 

0.478 

65.224 

84,841 

17 

130.000 

172.672 

173.080 

114.682 

-36.938 

33.704 

34.950 

1.246 

58,398 

71,888 

18 

135.000 

168.244 

168.635 

117.585 

-32.779 

22.593 

23.806 

1.213 

51.050 

56.585 

19 

140.000 

163.506 

162.995 

120.063 

-28.874 

10.000 

8.206 

-1.794 

42.933 

37.080 

^uvg 

= 0,622 '' 

^rms 

= 0.768" 








^12 

= 20.000 

cm ^22 

= 26.175 

cm 

L31 = 35. 

000 cm 

L4J = 

10.000 . 

cm 


^ 5 a 

= 7.765 

cm Lsh 

= 35.000 

cm 

Lea ^ 12. 

226 cm 

^ 6 b ~ 

40.000 

cm 


Vc 

= -142.000" 7^ 

= -154.000" 








4n 



Table 4,3 Generated displacements for the Stephenson-I linkage in Example 


when (< 95 ^/) of Table 4.1 are used as input values. 


i 

^I2i 



{^22i)g 

{^31i)g 

{^4u)a 



Mu 


1 

50.000 

-170.510 

-170.510 

55.485 

-75.948 

10,000 

9.719 

-0.281 

45.996 

85.667 

2 

55.000 

-169.056 

-169,056 

58.998 

-73.751 

22.593 

22.451 

-0.142 

48.054 

83.798 

3 

60.000 

-168.076 

-168.076 

62.622 

-72.154 

33.704 

33.493 

-0.211 

50.698 

74.353 

4 

65.000 

-167.514 

-167.514 

66.338 

-70.910 

43.333 

43.119 

-0.214 

53.852 

65.971 

5 

70.000 

-167.315 

-167.315 

70.124 

-69.830 

51.481 

51.444 

-0.037 

57.439 

58.726 

6 

75.000 

-167.434 

-167.434 

73.962 

-68.748 

58.148 

58.451 

0.303 

61.396 

52.802 

7 

80.000 

-167.836 

-167.836 

77.834 

-67.511 

63.333 

64.010 

0.677 

65.670 

48.480 

8 

85.000 

-168.496 

-168.496 

81.723 

-65.981 

S7.031 

67.929 

0.892 

70.219 

46.090 

9 

90,-000 

-169.397 

-169.397 

85.615 

-64.057 

69.259 

70.054 

0.795 

75.012 

45.889 

10 

95.000 

-170.534 

-170.534 

89.493 

-61.699 

70.000 

70.374 

0.374 

80.027 

47.926 

11 

100.000 

-171.909 

-171.909 

93.344 

-58.932 

69.259 

69.042 

-0.217 

85.253 

52.026 

12 

105.000 

-173.536 

-173.536 

97.151 

-55.817 

61.031 

66.289 

-0.748 

89.313 

57.894 

13 

110.000 

-175.438 

-175.437 

100.895 

-52.422 

63.333 

62.325 

-1.008 

83.668 

65.253 

14 

115.000 

-177.653 

-177.653 

104.553 

-48.803 

58.148 

57.271 

-0.877 

11 .194 

73.925 

15 

120,000 

179.761 

179.761 

108.099 

-45.000 

51.481 

51.138 

-0.343 

71.662 

83.862 

16 

125.000 

176.716 

176.716 

111.493 

-41.039 

43.333 

43.802 

0.469 

65.223 

84.841 

17 

130.000 

173.080 

173.080 

114.683 

-36.947 

33.704 

34.934 

1.230 

58.396 

71.881 

18 

135.000 

168.635 

168.635 

117.587 

-32.789 

22.593 

23.781 

1.188 

51.048 

56.570 

19 

^avg 

140.000 

= 0.624'' 

162.995 162.995 

0.773® 

120.064 

-28.889 

10.000 

8.158 

-1.842 

42.930 

37.048 

^12 

= 20.000 

cru '^22 

= 26.174 

cm 

L3i= 35.1 

001 cm 

^41 = 

10.000 cm 


Vc 

= 7.764 cm Ljj 

= -142.000® 

= 35.001 cm 

= -154.000 ® 

^6.= 12. 

225 cm 

^65 = 

40.000 cm 



Table 4.4 Desired and generated displacements for the Stephenson~I 
in Example 2 for an output oscillation angle of more than 180" . 

linkage 

i 

dm 

{dsai)a 

(^5m)g 

{d22i)g 

(^3w)g 




Mu 

Mli 

1 

0.000 

45.000 

43.356 

55.950 

-66.338 

-5.000 

-11.170 

-6.170 

12.594 

55.168 

2 

30.000 

15.000 

15.525 

32.490 

-87.674 

-60.000 

-50.346 

9.654 

16-965 

37.328 

3 

60.000 

16.000 

15,715 

56.777 

-69.493 

12.000 

8.597 

-3.403 

41.063 

78,090 

4 

90.000 

18.000 

17.413 

82.407 

-57.903 

55.000 

53.041 

-1.959 

64.994 

69.056 

5 

120.000 

20.000 

20,053 

107.104 

-48.890 

90.000 

90.374 

0.374 

87.051 

40.736 

6 

150.000 

24.000 

24.236 

129.328 

-40.525 

118.000 

120.113 

2.113 

74.908 

19.362 

7 

180.000 

31.000 

31.214 

146.629 

-30.358 

130.000 

133.307 

3.307 

64.586 

16,335 

8 

-150.000 

42.000 

42.154 

156.341 

-20.849 

140.000 

134.205 

-5.795 

65.814 

24.947 

9 

-120.000 

56.000 

56.210 

158.169 

-14.847 

130.000 

130.580 

0.580 

78.041 

34.573 

10 

-90.000 

71.000 

70.474 

153.347 

-12,233 

118.000 

120.585 

2.585 

82.874 

47.182 

11 

-60.000 

82.000 

80.676 

140.984 

-14,929 

100.000 

102.604 

2.604 

60.308 

62.466 

12 

-30.000 

75.000 

78.220 

114-402 

-28.327 

70,000 

70.350 

0.350 

36.182 

81.323 

^avg 

= 3.241" 


^rms ■“ ^ ■ 

,181" 








= 4.500 cm 

L 22 = 5 . 

600 cm 

^31= 8. 

.037 cm 

^41 = 

2.364 1 

cm 


^5a 

= 9.267 cm 

II 

00 

292 cm 

00 

ii 

.284 cm 

II 

10.000 

cm 


Ic 

= 12 . 000 ® 


r, = i3. 

000 ® 
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TaDie 4,0 Generauea aispiacemenus lor cxiti ^ 


when (^sai) of Table 4 . 4 are used as input values . 


i 

hli 

i^Sai 


(^22i)g 

{'931/)^ 




j^ii 

ij-ii 

1 

0.000 

43.356 

43.356 

55.950 

-64.493 

-5.000 

-10.421 

-5,421 

12,594 

54.072 

2 

30.000 

15.525 

15.525 

32.490 

-85.426 

-60.000 

-47.716 

12.284 

16.965 

37.710 

3 

60.000 

15.715 

15.715 

56.777 

-67.387 

12.000 

9.928 

-2.072 

41.063 

77.316 

4 

90.000 

17.413 

17.413 

82.407 

-55.964 

55.000 

53.600 

-1.400 

64.993 

70.436 

5 

120.000 

20.053 

20.053 

107.104 

-47.096 

90.000 

90.246 

0.246 

87.051 

42.658 

6 

150.000 

24.236 

24.236 

129.328 

-38.793 

118.000 

119.283 

1.283 

74.908 

21.924 

7 

180.000 

31.214 

31.214 

146.629 

-28.799 

130.000 

132.883 

2.883 

64.586 

18,318 

8 

-150.000 

42.154 

42,154 

156.341 

-19.244 

140.000 

134.022 

-5.978 

65.814 

26.734 

9 

-120.000 

56.210 

56.210 

158.169 

-12.971 

130.000 

129.894 

-0.106 

78.041 

37.136 

10 

-90.000 

70.474 

70.474 

153.347 

-10.167 

118.000 

119.478 

1.478 

82.874 

50.355 

11 

-60.000 

80.676 

80.676 

140.983 

-12.851 

100.000 

101.325 

1.325 

60.308 

65.824 

12 

-30.000 

78.220 

78.220 

114.401 

-26.397 

70.000 

69.285 

-0.715 

36.182 

84.318 

^avg 

= 2.933 ^ 

e 

4.442'’ 







hi 

= 4.500 cm hi ^ 5.600 cm 

hi = 

.984 cm 

hl = 

2.476 cm 


ha 

= 9.267 cm ^sb ^ 8,118 cm 

^6a = ^ 

.284 cm 

h,- 

10.000 cm 


nc 

= 12.000'’ 


/a = 13.000'’ 








Table 4.6 Desired and generated co-ordinates of the path for the Stephensoi 
linkage in Example 3 . 


i 

(p.l 

(.cm) 

U). 

(cm) 

(■) 

fe). 

t) 

««• 

(cm) 

eyi 

(cm) 

1 

75.000 

74.769 

12.000 

10.593 

0.134 

-1.849 

2 

73.000 

72.507 

5,000 

2.126 

-0.265 

-3.673 

3 

67.000 

68.336 

18.000 

18.687 

1.012 

1.191 

4 

67.000 

67.155 

49,000 

50.781 

-1,495 

1.462 

5 

71.000 

70.935 

82.000 

83.060 

-1.310 

0.106 

6 

78.000 

77.702 

106.000 

107.256 

-1.549 

-0.773 

7 

85.000 

84.652 

126.000 

125.062 

1.332 

0.524 

8 

92.000 

90.042 

139.000 

138.648 

1.842 

-0,868 

9 

94.000 

93.107 

149.000 

148.870 

0.874 

-0.279 

10 

94.000 

93.869 

156.000 

155,863 

0.212 

0.152 

11 

94.000 

92.685 

161.000 

159.718 

1.940 

1.525 

12 

88.000 

89.853 

159.000 

160.548 

-2.569 

-1.614 

13 

84.000 

85.791 

157.000 

158.130 

-2.295 

-0.864 

14 

78.000 

81.234 

152.000 

151.802 

-2.723 

1.766 

15 

78.000 

77.163 

141.000 

140.652 

0.947 

-0.163 

16 

77.000 

74.437 

126.000 

123.565 

4.104 

-0,270 

17 

72.000 

73.109 

96.000 

97.386 

-1.873 

0.896 

18 

73.000 

73.126 

51.000 

53.771 

-2.722 

2.256 

^avg 

= 2.746 

cm 


^mts ” 2 • 

382 cm 
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Table 4.7 Generated displacements for the Stephenson-I linkage in 


Example 3 for path generation. 


i 

^I2i 


{^Sai)g 


{^ni)g 


Mu 

Mu 

1 

0.000 

39.500 

39.502 

66.147 

- 78.123 

- 36.762 

26.645 

41.361 

2 

20.000 

19.600 

19.744 

57.338 

- 85.532 

- 47.626 

37.594 

37.906 

3 

40.000 

10.300 

10.328 

62.440 

- 64.820 

- 32.383 

52.113 

32.436 

4 

60.000 

5.600 

5.611 

73.242 

- 30.959 

2.418 

67.631 

33.377 

5 

80.000 

3.300 

3.266 

86.380 

- 2.876 

40.165 

83.113 

43.041 

6 

100.000 

2.500 

2.487 

100.330 

13.458 

70.960 

82,158 

57.502 

7 

120.000 

3.200 

3.155 

114.105 

22.525 

95.840 

69.049 

73.314 

8 

140.000 

5.600 

5.620 

126.737 

28.075 

116.042 

58.884 

87.967 

9 

160.000 

10.600 

10.595 

137.099 

32.458 

131.034 

53.495 

81.425 

10 

180.000 

18.700 

18.662 

144.287 

37.257 

139.502 

54.376 

77.756 

11 

- 160,000 

29.500 

29.547 

148.299 

42.744 

141.343 

61.248 

81.401 

12 

- 140.000 

42.200 

42.230 

149.855 

A1 ,169 

137.937 

72.376 

89,831 

13 

- 120.000 

55.600 

55.618 

149.588 

50.692 

130.530 

86.030 

79.838 

14 

- 100.000 

68.700 

68,723 

147.686 

49.878 

119.539 

78.964 

69.661 

15 

- 80.000 

80.400 

80.364 

143.763 

43.781 

104.763 

63.399 

60.982 

16 

- 60.000 

88.500 

88.456 

136.474 

30-957 

85.263 

48.018 

54.307 

17 

- 40.000 

88.300 

88.246 

122.341 

8.523 

57.006 

34.094 

48.483 

18 

- 20.000 

70.500 

70.429 

95.440 

- 32.623 

9.709 

25,011 

42.332 

^12 = 

60.000 cm 

^22 

= 85.674 cm 

^ 31 = 58 . 

.333 cm 


51.704 cm 

^5a- 

89.685 cm 

Lst 

= 78.949 cm 

it 

.459 cm 


21.863 cm 

Rj^ = 

5.316 cm 

a 

= 40.533 cm 





rc = 

- 20 . 000 *’ 


1 = 180 . 

000 " 

ap = - 103 . 000 " 

as = 

0.000 " 


Table 4.8 Desired and generated co-ordinates of the path for the Watt-H 
linkage in Example 4 for motion of the garage-door. 



U). 

tel tel 

^xi ^yi 

(cm) 

(cm) 

(') (') 

{cm) {cm) 


1 

41 . 

.000 

42 . 

684 

108 . 

000 

2 

40 . 

000 

39 . 

911 

102 . 

838 

3 

39 , 

,100 

37 . 

910 

94 , 

,100 

4 

39 , 

,100 

37 . 

423 

85 , 

,900 

5 

40 , 

.000 

38 . 

354 

77 , 

.162 

6 

41 . 

.000 

40 . 

085 

72 . 

.000 

7 

43 , 

.000 

42 . 

,466 

65 , 

.091 

8 

45 , 

.000 

46 , 

.339 

60 . 

.073 

9 

47 . 

.000 

49 , 

,389 

56 , 

.077 

10 

51 

.000 

54 . 

.182 

49 

.881 

11 

54 

.000 

59 . 

.038 

46 , 

.238 

12 

57 

.000 

62 , 

.538 

43 

.173 

13 

62 

.000 

63 , 

.146 

38 

.979 

14 

66 

.000 

64 , 

.722 

36 

.221 

15 

70 

.000 

66 . 

.225 

33 

.858 

16 

73 

.000 

69 , 

.441 

32 

.293 

^avg 

= 

3.075 

cm 





107 . 

718 

- 0 . 

320 

1 . 

666 

102 . 

190 

0 . 

460 

0 . 

Oil 

94 . 

,516 

- 0 . 

189 

- 1 . 

208 

85 , 

,881 

- 0 . 

107 

- 1 , 

,674 

78 . 

.039 

- 0 . 

,939 

- 1 . 

.479 

71 . 

.707 

- 0 . 

.088 

- 0 , 

.934 

66 , 

.019 

- 0 . 

,851 

- 0 , 

.199 

59 , 

.357 

1 . 

.168 

0 , 

.868 

55 , 

.137 

2 , 

.002 

1 . 

.525 

49 , 

.366 

2 

.421 

2 . 

.117 

44 

.012 

5 

.110 

2 , 

.020 

40 

.248 

6 

.163 

1 . 

.406 

39 

.595 

0 

.461 

1 , 

.246 

37 

.895 

-2 

.171 

0 

.754 

36 

.266 

-4 

.733 

0 

.175 

32 

.772 

-3 

.321 

-1 

.412 


= 2.995 cm 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


('’^“)« fe-)- (‘’«^). fell), («,„) 


-114.455 

-111.763 

-108.772 

-105.923 

-103.549 

-101.680 

- 100.000 

-98.018 

-96.761 

-95.061 

-93.516 

-92.439 

-92.251 

-91.757 

-91.270 

-90.124 


8.092 
7.761 

6.093 
2.823 

-1.575 

-6.421 

-12.063 

-20.601 

-27.213 

-37.731 

-48.764 

-57.039 

-58.505 

-62.352 

-66.073 

-74.115 


8.095 

7.763 

6.094 

2.823 

-1.576 

-6.422 

-12.064 

-20.601 

-27.213 

-37.732 

-48.764 

-57.038 

-58.505 

-62.352 

-66.073 

-74.117 


66.970 

59.951 

51.045 

41.041 

31.133 

22.015 

12.583 

-0.361 

-9.717 

-23.853 

-38.013 

-48.294 

-50.087 

-54.747 

-59.183 

-68.405 


£_ 

67.794 

60.770 

51.876 

41.909 

32.058 

23.007 

13.659 

0.857 

-8.377 

-22.289 

-36.170 

-46.210 

-47.958 

-52.496 

-56.815 

-65.819 




g 


e. 


Mii 


Mil 


9.000 

6.000 
4.000 
0.000 

-4.000 

-9.000 

-14.000 

- 22.000 

-28.000 

-38.000 

-48.000 

-56.000 

-61.000 

-67.000 

-73.000 

-81.000 


7.108 

6.642 

4.815 

1.376 

-3.188 

-8.187 

-13.992 

-22.768 

-29.568 

“40,406 

-51.813 

-60.401 

-61.925 

-65.925 

-69.794 

-78.111 


-1.892 

0.642 

0,815 

1.375 

0.812 

0.813 

0.008 

-0.768 

-1.568 

-2.406 

-3.813 

-4.401 

-0.925 

1.075 

3.206 

2.889 



= 

mts 

2.110® 



Lji = 44.933 cm 

II 

37.614 cm 

■^^31 — 26.177 cm 

L^l = 22 . 574 cm 

L^a~ 22.108cm 

Ls,= 

0.426 cm 

“ 15.000 cm 

= 41 . 446 cm 

= 46.976 cm 

II 

64.681 cm 

a= 76.221 cm 

i/ = 9.000 ® 

y,.= -20.000® 

Vc = 

0.000 

= 0.000 "" 

Oa = 44.000® 


57.450 

60.474 

65.134 

71.254 

78.027 

84.742 

87,937 

77.416 

69.548 

57.329 

44.752 

35.401 

33.746 

29.405 

25.198 

16.007 


60.687 

54.127 

47.061 

40.534 

35.246 

31.194 

27.651 

23.625 

21,192 

18.117 

15.643 

14.191 

13.967 

13.429 

12.979 

12.292 


5= 351.000' 


Table 4.10 Generated co-ordinates of the path for the Watt-II linkage in 
Example 4 for motion of the garage-door, with Lsi=(icm.. 


'■ M, (pi), fe). 



(cm) (cm) (®) 

("I (cm) 

(cm) 


1 

41. 

.000 

41. 

.430 

108, 

.000 

108. 

.319 

-0, 

.352 

0. 

.337 

2 

40. 

.000 

38. 

.621 

102. 

.838 

102. 

.669 

0. 

.418 

-1. 

.319 

3 

39, 

.100 

36. 

.593 

94. 

.100 

94, 

.813 

-0. 

.274 

-2 

.536 

4 

39. 

.100 

36. 

.097 

85. 

.900 

85. 

.990 

-0. 

.271 

-2 

.992 

5 

40. 

.000 

37. 

.032 

77, 

.162 

78. 

.029 

-1. 

.207 

-2 

.774 

6 

41. 

.000 

38. 

.773 

72, 

.000 

71, 

.652 

-0, 

.465 

-2 

.192 

7 

43. 

.000 

41. 

.168 

65, 

.091 

65. 

.972 

-1. 

.348 

-1 

.400 

8 

45. 

.000 

45. 

.066 

60, 

.073 

59, 

.374 

0. 

.508 

-0 

.219 

9 

47. 

.000 

48. 

.144 

56, 

.077 

55, 

.220 

1. 

.233 

0 

.543 

10 

51. 

.000 

53, 

.000 

49, 

.881 

49. 

.552 

1. 

.521 

1 

.333 

11 

54. 

.000 

57, 

.953 

46. 

.238 

44. 

.287 

4. 

.136 

1 

.466 

12 

57. 

.000 

61, 

.549 

43. 

.173 

40. 

.569 

5, 

.184 

1 

.030 

13 

62. 

.000 

62. 

.175 

38, 

.979 

39, 

.922 

-0, 

.514 

0 

.901 

14 

66. 

.000 

63. 

.804 

36. 

.221 

38, 

.236 

“3. 

.129 

0 

.489 

15 

70, 

.000 

65. 

.361 

33. 

.858 

36, 

.614 

-5, 

.666 

-0 

.017 

16 

73. 

.000 

68, 

.709 

32. 

.293 

33. 

,118 

-4, 

.161 

-1 

.460 


egyg= 3.212 cm 3.079 cm 
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Table 4.11 Generated displacements for the Watt-II linkage in Example 4 
for motion of the garage-door, with L 5 ^= 0 cm, 


”1 ^TuTasT 

2 -111.803 

3 -108.809 

4 -105,957 

5 -103.579 

6 -101.707 

7 -100.024 

8 -98.036 

9 -96.775 

10 -95.068 

11 -93.512 

12 -92.424 

13 -92.233 

14 -91.731 

15 -91.234 

16 -90.049 

1-853 " 




7.000 

8 

7.000 

7 

6.000 

6 

4.000 

2 

0.000 

-1 

-6.000 

-6 

-11.000 

-12 

-20.000 

-20 

-30.000 

-27 

-40.000 

-37 

-50.000 

-48 

-55.000 

-57 

-58.000 

-58 

-63.000 

-62 

-67.000 

-66 

-70.000 

-74 


^rms 



.014 66.970’ 
.686 59.951 
.022 51.045 
.756 41.041 
.639 31.133 
.483 22.015 
.123 12.583 
.661 -0.361 
.274 -9.717 
.798 -23.853 
.841 -38.013 
.131 -48.294 
.601 -50.087 
.460 -54.747 
.196 -59.183 
.305 -68.405 

2.174° 


(^ 31 /)^ 

~~ei .802 9“ 

60.781 6 

51.892 4 

41.931 0 

32.087 -4 

23.043 -9 

13.705 -14 

0.920 -22 

-8.298 -28 

-22.181 -38 
-36.024 -48 

-46.030 -56 

-47.772 -61 

-52.295 -67 

-56.599 -73 

-65.587 -81 


000 8.108 
000 7.643 

000 5.816 

000 2.374 

000 -2.194 

000 -7.200 

000 -13.015 

000 -21.811 
000 -28.631 

000 -39.508 

000 -50.971 

000 -59.612 

000 -61.147 

000 -65.177 

000 -69.078 

000 -77.478 




-0.892 

1.643 

1.816 

2.374 

1.806 

1.800 

0.985 

0.189 

-0.631 

-1.508 

-2.971 

-3.612 

-0.147 

1.823 

3.922 

3.522 


57.490 

60.511 

65.169 

71.287 

78.060 

84.776 

87.901 

77.375 

69.501 

57.270 
44.671 
35.293 
33.633 

29.271 
25.038 
15.743 


59.695 

53.138 

46.076 

39.557 

34.281 

30.243 

26.720 

22.731 

20.333 

17.327 

14.947 

13.582 

13.375 

12.883 

12.479 

11.892 


Lji = 44.933 cm 
^s« = 22 . 108crri 
R^= 46.976 cm 

- 20 . 000 ° 


L 21 = 37 - 614 cm 

11 

26.177 cm 

II 

22,574 cm 

L52, = 0.000 cm 

ii 

15.000 cm 

II 

41.446 cm 

Rb = 64.681 cm 


76.221 cm 


9.000 " 

Tj^= 0.000° 

«A = 

0.000° 

ag = 

44.000° 
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Chapter 5 


Optimization of Transmission Angle for 
Planar Four-Bar Linkages 

5.1 Introduction 

As stated in Chapter 1, the present chapter is an extension of the work done by Dubey [1] 
which covered the function generation and path generation problems for planar four-bar 
linkage. In this chapter, the solutions obtained in reference [1] have been optimized from 
the view' point of best transmission angle. Before going into the details of the optimization 
method, a brief note on transmission angle and its usefulness is given. 

5.2 Transmission Angle 

Before actually building a linkage, it is useful to have a measure of how well a linkage 
may perform while it is still in kinematic skeleton form on the drawing board or on the 
computer screen. Hartenberg and Denavit [7] gave a term "run" which formally means 
the effectiveness with which motion is imparted to the output link. It implies a smooth 
operation, in which a maximum force component is available to produce a torque or a 
force, whatever the case might be, in the output link. Though, a test on a kinematic 
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model checks out only an approximation of the static forces and tells nothing about the 

dynamic forces, however, some evaluation is better than none. IVansmission angle is taken 
as a measure of the run. 

Alt defined the transmission angle, ij,, as the smaller (acute) angle between the direc- 
tion of the velocity-difference vector cb of the floating link and direction of the absolute 
velocity vector of the output link both taken at the point of connection. The transmis- 
sion angle is illustrated in Figure 5.1. Since the velocity vectors are perpendicular to their 
respective links, the transmission angle is also given by the angle between the floating Hrik 
(iz) ^od the output link (^ 4 ) as shown in Figure 5 . 2 . Clearly the optimum value of fj. is 
90®; the recommended range is jj. >30®. 



Hence, it is always desirable to have the value of transmission angle more than 30°. 
The next section describes the optimization process as applied to the function generating 
linkages. 

5.3 Optimization for Function Generation 

The four-bar linkage for function generation is as shown in Figure 5.3. Link 2 is the input 
link and link 4 is the output link. The governing loop-closure equation is given by 


47 



c 


Figure 5.2. Transmission angle fi as the acute angle between link 3 and link 4. 


y 



Figure 5.3: Four-bar linkage at the F'* design position when generating a function. 
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(5.1) 


Il + L4-I)’2-^3 = 0 . 

Equation (5.1) after elimination of ^ 3 , at design position gives 

L\ + Ll-}-Ll-Ll- 2.Li.l2Cos(^2i) + 2 .Li.L^cos{ 9 ii) - 2.L2.L, cos(02i - ^4:) = 0 . (5.2) 
Dividing Equation (5.2) by 2.L2.D4 one gets 

DiHii + D2ii2i + Da/isi - 1 /^ = 0 (5.3) 

where, 

Di = 7^, jDs = Dz = 

1 £,4’ ^ i2’ 3 2.I,2.L4 > 

= cos (02i), /£2i = - cos (04i), ^3^ = -1, 

Vi = - cos (02j - ^4i) • 

The value of L\ is specified. 621 and On are obtained from the given function to be 
generated (i.e., x and y — f{x)). 

Referring to Figure (5.3), it is seen that 7 is the angle which ^4 makes with ~tz- R 
follows that, 

BD'^ = AB^ + AD^ - 2AB AD cos (02i) = BC^ + DC^ - 2. BC. DC cos ( 7 ^) 


or 


cos (7i) = 


BC^ + DC^ - AB^ - AD^ + 2 ABAD cos {621) 
2.BC.DC 


or 


cos (7,) = 
which gives (using Equation ( 5 . 3 )), 


L\a L\- Ll- L\ + 2.L2.L1 cos {92i) 
2.L3.L4 


cos (7i) 


2.L4 — 2.L2.L4.D3 + 2.L2 -Li cos (02i) 
2.LZ.L4 
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Sinrr /. !S liir afutr angle between and it follows that, 


30° < 7i < (180° - 30°) 


or 


30° < 7t < 150° 


which can he written as 


7i - 30° > 0 


and 

-7i + 150°>0. 

It should be noted that L3 is given by yjL\ + LI + LI- 2 .L 2 .L^.Dz (from Equation 
( 5 . 3 )) and for L3 to have a real magnitude, the terms under the square root should give a 
positive value. Taking, = 1,2 4. — 2.L2.T4.D3, the condition is written as. 


x>o. 

The optimization problem, is formulated as follows: 

Minimize 

N 

f D2, D 3 ) = Y^{DifJ,u + D2ji2i + Dzli-Zi — ^i) ( 5 - 4 ) 

1=1 

subject to 

D2 1^0, 

li - 30 ° > 0, 

-7i + 150 ° > 0, 





x> 0 , 


wluH i Hpioscnts the number of the design positions and all other variables have their 
designated meanings. The first two constraints are necessary else (with reference to Equa- 
tion (5.3)) L 2 and £4 will become indefinite. If the Grashofian condition is to be satisfied 
and £2 is to be the crank, then the following constraints are also included, 


£4 — £2 > 0, 


£3 — £2 > 0, 

£4 — £2 > 0, 

£1 -f £2 + £3 -h £4 — Z (£2 + Ljnax .) > 0 . 

The objective function is a linear function, but the constraints, particularly, the con- 
straint on the transmission angle, turn out to be highly nonlinear, making the optimiza- 
tion problem a nonlinear problem. Two approaches have been used for optimization, the 
Cauchy’s steepest descent method and the Random search method. 

Both the optimization methods need a start value i.e., some initial values of Di, D 2 , D 3 . 
The initial values are provided by solving the objective function without considering the 
constraints, using LSM. The search region for Di,D 2 ,Dz, when the constraints are in- 
cluded, has many local optima. This implies that the solution for the constrained problem 
will tend to localize near the initial values, but, this is certainly a better approach than 
giving arbitrary initial values. 

Moreover, the limit on optimization of the transmission angle is not kept 30° at start- 
ing. This is explained as follows. Suppose that a linkage obtained from solving the un- 
constrained problem using LSM results in a minimum transmission angle of 17 . The 
constraint put on the transmission angle is then the nearest higher multiple of 5, which 
in this case is 20° . The reason for this stems from the fact that a higher t^alue of trans- 

,, SNTRM. uBRMT' 


niission as tlu ( oiisuaint, makes the optimization unstable. Unstable optimization 
means, random wandering of the search direction in case of the steepest descent method 
and insutficicnt region and insufficient number of random solutions in case of the Random 
search method. 

I he optimization is carried out till an optimum value of 30 ° or more is reached or up to 
the \alue , after which no improvement is possible. A limit is always set on the maximum 
number of iterations as well as the maximum number of random generations. 

The basic algorithms of the Cauchy’s steepest descent method and Random search 
method are available in reference [ 14 ]. These algorithms were then modified for the linkage 
synthesis application. In the steepest descent method, the Bracket Operator Penalty is 
used. The multi-variable constrained problem is then reduced to a single variable uncon- 
strained problem, which is solved by using the Bounding Phase Method and the Golden 
Section Search Method. In the random search method, random values of Di,D2,Ds are 
generated near the region of the initial values of Di, D2, D3. The search region is increased, 
discarding the previous search region, up to a particular limit in case the solution is not 
found. 

A point to be noted that, if the initial values of Di, D2, Dz (obtained by the LSM) gave 
good function generation results (i.e., very less value of error in function generation), then 
there is a great chance that the function generation error will increase after maximizing 
the transmission angle. 

A question arises that in cases where a crank was required for input, why the Grashofian 
condition was not incorporated during the LSM. The answer for this is that in the LSM, 
inclusion of constraints is not at all possible. If the solution is non- Grashofian, it is ne- 
glected and the LSM is again applied with changed initial values. But in the optimization 
of transmission angle, whether it is Cauchy’s steepest descent method or Random search 
method, the values of Di, D2, Dz and hence the link-lengths continuously change during 
the search. It is possible for the search to enter into the non-Grashofian zone. Therefore, 
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the Grashofian condition is taken as a constraint to penalize the objective function. It is 
possible that, within the limit on the maximum number of iterations, no solution emerges 
or, within the region of search, no solution exists. 

5.4 Optimization for Path Generation 

This section gives the optimization problem formulation only, without going into the de- 
tails. The four -bar linkage for path generation is shown in Figure 5.4. There are two 



Figure 5.4; Four-bar linkage at the i*'‘ design position when generating a path, 
loops, OAB PO and ABC FA from which the following equations are derived, 

Difiii + D2lMii + Dzfxzi - DiDzPi — I'i = 0 
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\vh('r(\ 


Dj /?,!, D-2 - 1^2, = 1^2 + - a^ A = DiD-2 = Ra.L2 

/iu = -~2.p, cos (a,i - 5i), 

Il2r = - 2 . Pi cos (^2i - ^i). M3i = 1. 

/;, = -2 cos {a A - 6'2i), Ui = -P^ 
and 


d" ^2^21 d~ ^zt^Zi ~ ^1^2pi 


z /,=0 


( 5 . 6 ) 


where, 

d; = . 4 = Rs , d; = Ri + LH Ro + Ri - , 

~ 2.Ra cos (q;^ o;£)) 2.L2. cos (^2i ~ ct£>). 

/l£2i = 2.R4 cos {aA - 6zi) + 2.L2 cos (02i - 03i), /^si = 1 , 

Pi = 2 cos (^3i - ao), 

v'i = -R4.R2COS (aA - 02i) . 

aA and a^ are assumed and 621, l^i (Pi with 6i) are specified. Using LSM, the uncon- 
strained problem is solved to give the values of Di,D2,Dz and D[, D'2,D'^ and hence the 
link-lengths. 

For the constrained problem, the constraints are just the same as that for the func- 
tion generation problem. The objective function constitutes only Equation ( 5 . 6 ) and not 
Equation ( 5 . 5 ). An objective function constituting both the equations could be taken and 
represented as follows: 

Minimize 

'tf{DuD^,D,) + Z/(D[,D'^,D',), 

i=l i=l 

but then it would be a problem of six variables. Moreover there would be unnecessary 
effects of E/ and E/' on each other, making the optimization process very complex. 
Further, taking Equation ( 5 . 5 ) as objective function will change the length of L2 which is 
obtained minimizing the error in the loop OABPO involving fiS. Any change in L2 can 
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■" P"’>‘ Therefore, only Equation (5.6) is talren as the 

objective function and one can write 

Minimize 



Ihe initial values of Dj, D^. are obtained from the unconstrained LSM solution. Both 
the methods of optimization, the Cauchy’s method and Random search method are used 
for path generation also. 

5.5 Numerical Examples 

Two examples are given which show the application of the optimization of transmission 
angle on the four-bar linkage. 

Example 1: A crank-rocker example was illustrated in Dubey [1] (p. 50), to show 
the application of least square method to function generation of four-bar linkage. The 
synthesized linkage had a very poor minimum transmission angle fimin = 17.685°.^ The 
values of the error in the generated function were given by eavg = 5.250'’ and Cms = 6.914°. 

The same data is taken as the input for the present program. The value of Li is taken 
as 10 cm. Table 5.1 shows the desired and generated values of the displacements of the 
resulting linkage. The Cauchy’s steepest descent method and Random search method are 
applied separately to maximize fXmin- Table 5.2(a) gives the values of the displacements 
generated and the dimensions obtained when the linkage is optimized for pmin = 20°, using 
Cauchy’s steepest descent method. The value of prnin obtained is 19.995°. It is seen that 
there is change in errors generated at each design position. In Table 5.2(b) the linkage 
is optimized for 25° and in Table 5.2(c) for 30°. The final values of error obtained are 

^Due to different program codings and different numerical accuracy maintained therein, the values of 
angles generated by the present program and Dubey’s [1] program differ upto a maximum of 3°. Dubey 
[1] (Plot 3.4, p. 60) reported fimin = 14.83°. 
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(■av^ = “ "94" and = 9.008". The optimized value of is 30.000". Tables 5.3(a) and 
5.3(1)) give the values obtained using Random search method. No further improvement is 
obtained after = 25.245". The solutions obtained by both the methods differ from 
ea( li other. Moreo\er. time taken by Random search method is more as compared to 
Cauch) s method. I he leason for this is the large number of search points generated in 
the Random .search method. 

It follows that both the methods improve the minimum transmission angle of the linkage 
at the same time increasing the error in the function generated. At this stage it can not 
be said which method is the better one. If Random search method fails to improve 
up to 30", then Cauchy’s method gives a large error of 18.599°, at the design position, 
for IJ.min = 30". 

Example 2: A four-bar linkage, shown in Figure 5.5, is used to generate a path whose 
coordinates are given in Table 5.4. Crank AB is the input link. The specified values are 
aA = 120", ao = 25" and V’ = -22". The errors in the synthesized linkage are given by 
&avg — 0.025 cm and erms = 0.052 cm. The transmission angle has a minimum \alue of 
Pmin = 23.249". 

Tables 5.5(a) and 5.5(b) give the values of the displacements and firnin obtained using 
Cauchy’s steepest descent method. The final values obtained are Cavg = 0.052 cm, e^ms = 
0.138 cm and /Zmin = 30.000". Tables 5.6(a) and 5.6(b) give the values obtained using 
Random search method. The values obtained are eavg = 0.046 cm, Crms = 0.118 cm and 
l^min - 30.147". Both the methods successfully optimized the minimum transmission angle 
but with increased error in the path generated. 

From Table 5.5(a) it is seen that application of Cauchy’s method for pram = 25° results 
in decrease in values of errors to Cavg = 0.009 cm and Crms = 0.009 cm. 
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Figure 5.5: Four-bar linkage synthesized in Example 2 for path generation. 

5.6 Conclusion 

The optimization of the transmission angle leads to a linkage with acceptable aptness of 
run. 

An ideal case would have been a high minimum transmission angle resulting from 
the unconstrained equation itself. But due to the initial assumptions and the various 
mathematical input-output relationships, this condition is attained in a very few cases 
only. The better option is to have more than one solution using LSM (by varying the set of 
initial assumptions) and subjecting them to optimization of transmission angle. The best 
solution obtained thereafter is the optimum solution. 

The satisfactory performance of the optimization methods largely depend on the math- 
ematical nature of the generated function and the arbitrary choices - starting angles and 
positions of hinges etc. When, even after application of both the methods, a poor trans- 
mission angle results, the only remedy is to redesign on different arbitrary choices. 

The case of motion generation, though not considered in the present work, is just similar 
to path generation. The only difference is that in the loop OABPO of Figure 5.4 instead 
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of (h- f’liniinatc'ci. 


Table 5.1 Desired and generated displacements for the four-bar linkage 
in Example 1 icr function generation. 


(^ 4 .-), e, 


1 

15 . 

000 

13, 

,852 

30. 

,000 

31, 

.538 

2 

60. 

000 

12, 

.990 

50, 

0 

0 

0 

47, 

.814 

3 

i05. 

000 

18, 

.650 

70, 

,000 

71, 

.334 

4 

150. 

ooc 

27. 

.153 

90. 

,000 

91, 

.405 

5 

-165. 

^000 

37, 

.854 

110, 

,000 

103. 

.891 

6 

-120. 

000 

48, 

.393 

100, 

,000 

105. 

.848 

7 

-75. 

,000 

52, 

.303 

85, 

.000 

93, 

.480 

8 

-30. 

,000 

37, 

.486 

75. 

.000 

59, 

.897 

e =5 

.250 ' 

* 


^nns 

= 6.914® 



= 1. 

235 



D,= 

3.080 



= 10 

.000 

cm 


Lj = 

3.247 cm 

L4 

“ 8. 

095 cm 


Mmn ■ 

= 17 

.685' 

> 


1. 

,538 

17, 

.685 

-2. 

.186 

34, 

.824 

1. 

,334 

52. 

.684 

1. 

.405 

64. 

.252 

-6, 

.109 

66. 

.037 

5. 

.848 

57, 

.454 

8, 

.480 

41. 

.178 

-15. 

.103 

22, 

.411 


I>3= -0.473 
1 - 3 = 14.176 cm 


Table 5.2(a) Generated displacements for the four-bar linkage in Example 1 
obtained using Cauchy's steepest descent method for //^ = 20® . 


i 

^2i 






1 

15.000 

17.539 

30.000 

37,534 

7.534 

19.995 

2 

60.000 

15.213 

50.000 

50.309 

0.309 

35.096 

3 

105.000 

20.314 

70.000 

71.979 

1,977 

51,665 

4 

150.000 

28.559 

90.000 

91.051 

1.051 

62.492 

5 

-165.000 

38,988 

110.000 

103.150 

-6,849 

64.163 

6 

-120,000 

49.187 

100.000 

105.314 

5.314 

56.127 

7 

-75.000 

52.964 

85.000 

93.914 

8.914 

40.950 

8 

-30.000 

39.449 

75.000 

63.432 

-11.568 

23.983 

5-440 

0 


.617'' 




=: 1.189 


D,= 3. 

188 

I >3 = 

-0.456 


= 10.000 

cm 

00 

li 

136 cm 

il 

14.302 cm 


= 8.407 

cm 

19-995® 

time- 

1.480 s 
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Table 5.2 ib) Generated displacements for the four-bar linkage in Example 
obtained using Cauchy's steepest descent method for u ■ =2^® 

“min 


i 


(«>.), 

(««). 

M. 

e.- 

Mi 

I 

15.000 

18.762 

30.000 

43.759 

13.759 

24,997 

<c 

60 . 000 

15.424 

50.000 

54.476 

A. 4.-76 

39,052 


105.000 

19.104 

70.000 

74.683 

4,683 

55.580 

4 

150.000 

26.055 

90.000 

92.717 

2.717 

66.661 


-165.000 

35.223 

110.000 

103.610 

-6.390 

68.387 

€ 

-120.000 

44.178 

100.000 

104.304 

4.304 

60.126 

7 

-75.000 

47.166 

85.000 

91.979 

6.979 

44.813 

t 

-30.000 

36,002 

75.000 

64.539 

-10.461 

28.538 

€ 

= 6.72l‘ 


^rms ~ ^ ‘ 

SSl" 



D, 

= 1.389 


D2= 3.744 

Z>3= H 

0.485 

L,= 

= 10.000 

cm 

1/2 = 2 . 671 cm 

1.3= 13.326 cm 

L 4 : 

= 7.197 

cm 

Mjmn 

.997® 

time = ; 

1.540 s 


Table 5.2(C) Generated displacements for the four -bar linkage in Example 
obtained using Cauchy's steepest descent method for /i^=30‘’. 


‘ 


1 

15.000 

18.599 

30.000 

48.599 

18.599 

30.000 

2 

60.000 

14.982 

50.000 

58.194 

8.194 

43,212 

3 

105.000 

17.555 

70.000 

77.199 

7.199 

59.644 

4 

150.000 

23.330 

90.000 

94.270 

4.270 

70.940 

5 

-165,000 

31.277 

110.000 

103.990 

-6.010 

72.714 

6 

-120.000 

39.021 

100.000 

103.278 

3.278 

64.256 

7 

-75.000 

41.343 

85.000 

90.211 

5.211 

48.868 

8 

-30.000 

32.190 

75.000 

65.406 

-9,594 

33.216 

■7-794 

0 

9 - 008 " 



Dt 

= 1.640 


4.450 

1>3 = 

-0.526 


= 10.000 

cm 

L2 = 2.247 cm 


12.515 cm 


= 6.097 

cm 

/<nvu.= 30 

.000® 

time- 

2.370 s 
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Table 5.3(a) Oenarated displacements for the four-bar linlcage in Example 

obtsiriGd. usinQ" Rsiidom sssircii itiGtlioci foir u. • =20^ 

“nun 



1 15.000 15.540 30.000 37.562 7.562 22.021 

2 60.000 12.938 50.000 52.324 2.324 39.386 

3 105.000 17.241 70.000 75.786 5.786 58.545 

4 150.000 24.762 90.000 96.043 6.043 71.281 

5 -165.000 34.917 110.000 108.187 -1.813 73.269 

6 -120.000 45.362 100.000 109.129 9.129 63.767 

7 -75.000 49.636 85.000 95.761 10.761 46.125 

8 -30.000 36.774 75.000 63.395 -11.605 26.621 

6.878 " e^= 7.662" 

Dj = 1.415 ©2= 3.259 D^= -0.294 

Lj = 10.000 cm L 2 = 3.068 cm Lj = 13.119 cm 


^-4 = 7 . 068 cm 


22 . 021 ' 


tme= 23.170 s 


Table 5.3(b) Generated displacements for the four -bar linkage in Example 
obtained using Random search method for //n„n = 2S®. 


i 

dii 




K 

). 




Mi 

1 

15. 

000 

16. 

246 

30, 

.000 

41. 

491 

11. 

491 

25 

.245 

2 

60. 

000 

12. 

715 

50. 

.000 

55. 

561 

5. 

561 

42 

.846 

3 

105. 

,000 

16. 

106 

70, 

.000 

79. 

094 

9. 

094 

62 

.988 

4 

150. 

,000 

22. 

915 

90, 

.000 

99. 

564 

9. 

564 

76 

.649 

5 

-165, 

,000 

32. 

662 

no, 

.000 

111. 

465 

1. 

465 

78 

,803 

6 

-120. 

.000 

43. 

030 

100 

.000 

111. 

592 

11. 

592 

68 

.562 

7 

-75. 

.000 

47. 

538 

85 

.000 

97. 

409 

12. 

409 

49 

.872 

8 

-30. 

.000 

35. 

925 

75 

.000 

65. 

726 

-9. 

274 

29 

.800 

^avg 

, = 8 

.806 

o 


^rms ’ 

= 9. 

445'' 







= 1. 

560 



II 

3.405 


D 

3 = 

-0.148 

L, 

= 10 

,000 

cm 


.r 

II 

2.937 cm 


L 

3 ” 

12,462 cm 

L, 

= 6 

.410 

cm 


Anin 

= 25 

.245" 


time^ 

13, 

790 s 


No further improvement in within the specified limits 

Of the optimization zone. 
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a<=xxcxcn.ea aispiacements for the four-bar linkage i 
Example 2 for path generation. 


I 

(cm) 

(cm) 

(4). 

(”) 

« 

««• 

(cm) 

(cm) 

('’) 

(^l), 

(•) 

(') 

A 

(•) 

1 

3.779 

5.787 

33.790 

33.752 

0.008 

0.002 

-36.000 

40.313 

64.845 

36.018 

2 

3.994 

5.887 

25.670 

23.612 

-0.039 

-0.173 

-6.000 

6.718 

33.584 

23.249 

3 

3.867 

3.863 

34.050 

34.145 

-0.007 

0.003 

24.000 

-5.294 

38.354 

29.448 

4 

3.694 

3.692 

49.830 

50.020 

-0.011 

0.006 

54.000 

-7.896 

55.394 

47.337 

5 

3.560 

5.558 

67.970 

67.675 

0.016 

-0.009 

84.000 

-6.827 

74.944 

66.970 

6 

3.272 

3.274 

84.440 

84.529 

-0.005 

0.002 

114.000 

-2.942 

93.567 

84.943 

7 

2.969 

:.970 

97.860 

97.889 

-0.002 

0.001 

144.000 

4.301 

108.850 

81.354 

8 

2.588 

2.588 

105.370 

105.352 

0.001 

0.000 

174.000 

15.459 

118.966 

74.913 

9 

2.198 

2.198 

103.870 

103.876 

0.000 

0.000 - 

•156.000 

29.783 

123.516 

77.743 

10 

2.000 

1.999 

90.830 

90.898 

-0.002 

-0.001 - 

■126.000 

45.053 

122.987 

88.901 

11 

2.246 

2.243 

71.110 

71.157 

-0.003 

-0.002 

-96.000 

57.810 

116.787 

74.497 

12 

2.960 

2.962 

52.550 

52.538 

0.002 

0.001 

-66.000 

61.221 

100.700 

55.224 

^avg 

= 0.025 cm 


0. 

052 cm 






D'l- 

= -4.23' 

/ 

Z>2 = 

2.285 

D'i- 

= 25.272 






: 3.901 

cm 

1-2 = 

2.180 cm 

2-3 = 

: 4.237 cm L, 

1= 3. 

3 86 cm 


a = 

2.983 

cm 

a A = 

120.000® 

aj) = 

= 25.000® 

¥ 

= -22. 

000® 



= 2.969 

cm 

Rd ~ 

2.285 cm 

Mtmn 

= 23.249 

0 





Table 5.5(a) Generated displacements for the four -bar linkage in Exairple 
obtained using Cauchy's steepest descent method for ■ 


(Pi)a 

(cm) 




{cm) 



(cm) (cm) 


^2i 

(") 


(^x), 

(•) 


(««), 

(') 


Mi 

(•) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


avg 


3.779 

3.787 

33.790 

33.752 

0,008 

0.002 

-36.000 

43.673 

70.223 

37.657 

3.994 

4.008 

25.670 

25.852 

0.007 

0.017 

-6.000 

10.488 

38.630 

25.001 

3.867 

3.863 

34.050 

34.145 

-0.007 

0.003 

24,000 

-3.068 

41.622 

31.100 

3.694 

3.692 

49.830 

50.020 

-0.011 

0.006 

54.000 

-6.644 

58.031 

49.082 

3.560 

3.558 

67.970 

67,675 

0.016 - 

•0.009 

84.000 

-6.296 

77.495 

69.123 

3.272 

3.274 

84.440 

84.529 

-0,005 

0.002 

114.000 

-3.038 

96.299 

87.682 

2.969 

2 . 970 

97.860 

97.889 

-0.002 

0.001 

144.000 

3.707 

111.826 

77.988 

2.588 

2.588 

105.370 

105.352 

0.001 

0.000 

174.000 

14.653 

122.004 

71.171 

2.198 

2.198 

103.870 

103.876 

0,000 

0.000 - 

■156.000 

29.077 

126.394 

74.175 

2.000 

1.999 

90 . 830 

90.898 

-0.002 - 

■0.001 - 

■126.000 

44.638 

125.674 

85 . 905 

2.246 

2.243 

71.110 

71.157 

-0,003 - 

■0.002 

-96.000 

57.859 

119.473 

76.868 

2.960 

2.962 

52.550 

52.538 

0,002 

0.001 

-66.000 

62.223 

103.995 

57 . 105 

= 0.009 cm 


0- 

009 cm 






= -4.163 


D2 = 

2,376 

D3 = 

25.274 





= 3.961 

cm 

Lj = 

2.180 cm 


4.163 cm 

II 

U) 

357 cm 


2.983 

cm 

a A = 

120.000® 

ao = 

25.000® 


y/^ -22. 

000® 


= 2.969 

cm 

Rd = 

2.377 cm 

/^min ■■ 

= 25.001 

0 

time = 0. 

550 s 



fil 


ociieidcea dispiacements for the four-bar linkage in Exampl 

obtained using Cauchy's steepest descent method for . 


i 


{pih 

fe). 



Cyi 

On 

M. 


Mi 


(cm) 

(cm) 

(•) 

(•) 

(cm) 

(cm) 

(“) 

e) 

e) 

(i 

1 

3.779 

3.787 

33.790 

33.752 

0.008 

0.002 

-36-000 

50.600 

82.901 

42.698 

2 

3.994 

4.276 

25.670 

31.040 

0.064 

0,475 

-6.000 

19.344 

51.524 

30,000 

3 

3.867 

3.863 

34.050 

34.145 

-0.007 

0.003 

24.000 

2.209 

50.568 

36-021 

4 

3.694 

3.692 

49.830 

50.020 

-0.011 

0.006 

54.000 

-4.071 

65.548 

54.647 

5 

3.560 

3.558 

67,970 

67.675 

0.016 - 

■0.009 

84.000 

-5.912 

85.093 

76,283 

6 

3.272 

3.274 

84.440 

84.529 

-0.005 

0.002 

114.000 

-4.732 

104.832 

82.868 

7 

2.969 

2.970 

97,860 

97.889 

-0.002 

0.001 

144.000 

0-247 

121.494 

65,894 

8 

2.588 

2.588 

105.370 

105.352 

0.001 

0.000 

174-000 

10.566 

131.888 

57.350 

9 

2.198 

2,198 

103.870 

103.876 

0.000 

0.000 - 

156.000 

25.633 

135.385 

61.165 

10 

2.000 

1.999 

90.830 

90.898 

-0.002 - 

■0.001 - 

126.000 

42.372 

133.618 

75.399 

11 

2.246 

2.243 

71.110 

71.157 

-0.003 - 

■0.002 

-96.000 

57.060 

126.936 

84.871 

12 

2.960 

2.962 

52.550 

52.538 

0.002 

0.001 

-66.000 

63.836 

112.431 

63.213 

^avg 

= 0.052 cm 


= 0. 

138 cm 






d\ 

= -3.924 

Di = 

2.614 

I>3 = 

25.293 






=: 4.123 

cm 


2.180 cm 


3.924 cm 

L4= 3. 

242 cm 


a = 

= 2.983 

cm 

a A = 

120.000® 

CCj^ = 

25.000® 


li 

to 

10 

.000® 


Ra 

= 2.969 

cm 

Rd = 

2.615 cm 

/^min 

= 30.000 

0 

time = 1. 

.430 s 



Table 5.6(a) Generated displacements for the four-bar linkage in Exair^le 
obtai ned using Random search method for • 

{ Pt ), {<^*1 (<^4 ^ 

(cm) (cm) ('') ('') (cm) (cm) ('') ('') ('') ('') 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


3.779 

3.787 

3.994 

4.064 

3.867 

3.863 

3.694 

3.692 

3.560 

3.558 

3.272 

3.274 

2,969 

2.970 

2.588 

2.588 

2.198 

2.198 

2.000 

1.999 

2.246 

2.243 

2.960 

2.962 


33.790 

25.670 

34.050 

49.830 
67.970 
84.440 
97.860 

105.370 

103.870 

90.830 
71,110 
52.550 


33,752 

26.915 

34.145 

50.020 

67.675 

84.529 

97.889 

105.352 

103.876 

90.898 

71.157 

52.538 


0.008 

0.024 

-0.007 

- 0.011 

0.016 

-0.005 

- 0.002 

0.001 

0.000 

- 0.002 

-0.003 

0.002 


0.002 

0.110 

0.003 

0.006 

-0.009 

0.002 

0.001 

0.000 

0.000 

- 0.001 

- 0.002 

0.001 


-36.000 

- 6.000 

24.000 

54.000 

84.000 

114.000 

144.000 

174.000 
-156.000 
-126.000 

-96,000 

- 66.000 


44.989 

12.288 

-2.171 

-6.470 

-6.775 

-4.171 

2.033 

12.816 

27.473 

43.421 

57.081 

62.096 


73.264 

41.645 

43.865 

60.196 

79.950 

99.237 

115,205 

125.422 

129.420 

128.194 

121.608 

106.117 


39.324 

26.479 

32,643 

51,049 

71.816 

88.691 

73.359 

65,927 

69.215 

81.868 

79.914 

59.333 


^rmK == 0,033 cm 

ntxs 


^avg~ 0.018 cm 

J}\ =-4.074 
Lj = 4.009 cm 

a= 2.983 cm 
= 2.969 cm 


D2= 2.447 
L2 = 2 . 180 cm 

120.000'' 

Rjj^ “ 2 . 447 cm 


D’3 = 25.435 
L3 = 4.074 cm 

ao= 25.000® 

Amin= 26.479® 


L4 = 3.274 cm 

y/^ - 22 . 000 ® 

time -26.100 s 




uj.tipxacements for the four-bar linkage in Exampli 

2 obtained using Random search method for = 30" . 


i 


(. 



0. 

(. 

5,), 

e 

xi 

i 

->i 

Oli 






{cm) 

{cm) 

( 

■) 

1 

:•) 

{c 

m) 

( 

cm) 

( 

•) 

( 


(“ 

) 

( 

;) 

1 

3. 

.779 

3, 

.787 

33, 

790 

33. 

752 

0. 

008 

0. 

,002 

-36. 

000 

49. 

192 

81. 

,707 

43. 

,121 

2 

3. 

.994 

4. 

.236 

25, 

,670 

30. 

240 

0. 

059 

0. 

,403 

-6. 

,000 

17. 

969 

50. 

.334 

30. 

,147 

3 

3. 

.867 

3. 

,863 

34. 

,050 

34, 

,145 

-0. 

007 

0, 

,003 

24, 

.000 

1. 

,155 

50. 

.027 

36. 

.304 

4 

3 , 

.694 

3. 

.692 

49. 

,830 

50. 

,020 

-0. 

Oil 

0. 

.006 

54. 

.000 

-5. 

,025 

65, 

.547 

55. 

.312 

5 

3. 

.560 

3. 

.558 

67. 

.970 

67, 

,675 

0, 

,016 

-0, 

.009 

84. 

.000 

-6. 

,937 

85, 

.563 

77, 

.417 

6 

3. 

.272 

3. 

.274 

84. 

.440 

84, 

,529 

-0, 

005 

0. 

.002 

114. 

.000 

-5. 

,953 

105, 

.773 

81, 

.155 

7 

2, 

.969 

2, 

.970 

97. 

.860 

97, 

,889 

-0. 

,002 

0, 

.001 

144. 

.000 

-1. 

,196 

122, 

.826 

63. 

.502 

8 

2 

.588 

2. 

.588 

105, 

.370 

105. 

,352 

0, 

.001 

0, 

.000 

174, 

.000 

9. 

.109 

133, 

.228 

54, 

.479 

9 

2 

.198 

2, 

.198 

103, 

.870 

103. 

.876 

0. 

.000 

0 

.000 

-156, 

.000 

24. 

.379 

136 

.342 

58 

.524 

10 

2 

.000 

1, 

.999 

90, 

.830 

90, 

.898 

-0. 

.002 

-0 

.001 

-126, 

.000 

41. 

.300 

134 

.101 

73 

.420 

11 

2 

-246 

2 

.243 

71, 

.110 

71, 

.157 

-0. 

.003 

-0 

.002 

-96. 

.000 

56. 

.050 

126 

.940 

86 

.222 

12 

2 

,960 

2 

.962 

52 

.550 

52, 

.538 

0. 

.002 

0 

.001 

-66 

.000 

62. 

.723 

111 

.886 

64 

.055 

II 

0.046 ' 

cm 



^rms 

= 0. 

118 

cm 











Di 

= . 

-3.903 


D2 = 

2.618 



= 25.606 








Lv 


4.126 cm 

L 

'2 ^ 

2.180 

cm 


L3 = 

= 3 

.903 

cm 


L4 = 

= 3 

.170 

cm 



a 

= 

2.983 cm 

II 

120.000" 


«2) 

= 25.000" 



-22 

.000' 

0 



Ra 

= : 

2.969 

cm 

II 

Q 

2.618 cm 


/^min 

30.147® 


time = 11 

.090 

s 




Chapter 6 


Conclusions 


In the present thesis, a software package for the synthesis, analysis and animation of planar 
six-bar linkages has been prepared. The synthesis of the four-bar linkage has been extended 
to include the consideration of transmission angle. The main conclusions that can be drawn 
from this work are summarized below: 

1. The synthesis of the linkages have been carried out by using the Least Square Er- 
ror minimizing Optimization Method (LSM)and the Linear Superposition 
Technique. The LSM is an efficient method. There is no limitation on the number 
of the design points and the structural error is optimized over the whole range of the 
linkage motion. The absence of iterative processes is an advantage of the method 
over other optimization methods. The present thesis covers Stephenson-I linkage for 
the function and path generation, and Watt-II linkage for the motion generation. 
The other linkages, viz., Stephenson-II and III and Watt-I can also be used for the 
synthesis process. 

The LSM also has certain limitations. The major difficulty is finding a good set of 
values of the specified dimensions and the assumed values of the angular positions 
of the floating link (angle in all the three cases of the synthesis) at the design 


points. 



Towards this end, coupler curve atlases for the planar six-bar linkages designating the 
angular positions of coupler links at certain reference positions will be very helpful. 

2. The optimization of transmission angle for four-bar linkage gives a better linkage 
from the view point of motion and force transference. At the same time, error in the 
generated functions also increase. Both the methods, viz., Cauchy’s steepest descent 
method and Random search method should be applied and the better solution should 
be accepted. 

3. The dimensions of the linkages and the values of the error in the generated functions, 
known as structural error, were given up to the third place of decimal. An additional 
error, known as mechanical error, due to deflections of the links, play in the joints, 
and manufacturing tolerances, will inevitably occur in addition to structural error 
in any actual linkage. Hence, depending on the need and practicality, the level of 
accuracy in the dimensions should be taken. 

4. The analysis of linkages is based on the analytical approach where the loop closure 
equations of the linkages are used. Each equation gives rise to two scalar equations 
by taking the x and y components separately. The resulting equations are then used 
for displacement analysis. The time derivatives of the equations are used for velocity 
and acceleration analysis. The software provides many options for the analysis, which 
are given in Appendix A. 

5. The animation of the linkages provide a better understanding of the actual motion of 
the linkages. The stationary configurations and the displacement zone of the linkages 
can be easily seen from the animation. 

6. In Chapter 2, it was stated, and supported with examples, that six-bar linkages have 
many advantages over the four-bar linkages. Hunt [4] (pp. 224 - 34) cites certain 
examples of the six-bar linkages, where complex motion is generated by them, but 



with poor traiisniissioa angles. The planar six-bar linkages provide a fruitful field for 
further investigation. 


6.1 Suggestions for Future Work 

The following directions may be taken up to extend the work carried out in the present 
thesis: 

1. Weighted error function can be used in the synthesis processes, where the critical 
design points can have more weightage than the other design points. (The weightage 
error function for planar four-bar linkages was applied by Dubey [1].) 

2. The optimization of transmission angles for six-bar linkage can also be pursued. The 
problem becomes more complex than the four-bar linkage as there are two transmis- 
sion angles and more loops-closure equations involved. 

3. Tolerance Analysis can be included, which will give an indication of the magnitude 
of mechanical errors. In the software, the dimensions and angular positions of the 
linkages are given up to three places after the decimal point. This can be controlled 
up to the accuracy really possible. 

4. An atlas of the coupler curves of all the linkages with varying dimension ratios can 
be obtained, which can be treated as a guide for the designers. 

5. The dynamics of linkages should be included, so as to make the package more useful. 
The design of linkages depends to a great extent on the dynamics involved in the 
motion. 

6. The package should be ported to Windows environment, so as to have all the facili- 
ties and tools of a commercially available software package. This is being currently 
attempted as a part of the same DST project. 
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Appendix A 


The kinematic analysis of a linkage gives its performance characteristics. The analysis 
deals with the displacement, velocity and acceleration of the member links. 

The software package provides the following facilities for the kinematic analysis of a 
six-bar linkage: 

1. Displacement analysis: For the specified dimensions and input movement of a linkage, 
the movements of all other links are obtained. The ranges of angular displacements of all 
the links can also be seen. 

2. Velocity analysis: The angular velocity of any link (for a prescribed input) is obtained 
both in tabular and graph forms. 

3. Acceleration analysis: The package also gives the complete acceleration analysis of the 
linkage. 

The analysis is based on the analytical method of differentiating the loop closure equations 
of the linkage with respect to time. The details of the method can be referred to in [6]. 

Apart from these facilities, transmission angle curves can also be obtained. Provision 
for kinematic inversion is present by allowing the user to change the fixed link of the linkage. 
Coupler curves drawn by a specified point on the floating link can be obtained. This gives 
insight to the path behavior of the coupler point. 



